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ABSTRACT 


Sound  field  characteristics  were  studied  as  a  function  of 
increasing  intensity.  It  was  determined  that  one  effect  out¬ 
weighs  others.  This  is  the  distortion  generated  within  the 
siren  which  makes  it  impossible  to  generate  pure  tones  above  a 
certain  intensity.  At  high  intensity  the  waveform  of  the  sound 
assumes  a  sawtooth  shape  producing  periodic  shock  waves.  Con¬ 
versely,  it  is  concluded  that  if  the  distortion  is  small  at  the 
mouth  of  the  siren,  additional  distortion  in  the  subsequent 
propagation  will  not  be  significant.  Waveforms  observed  in  the 
AFFDL  facility  are  included  in  an  appendix. 

Interaction  of  sound  waves  from  different  sirens  will,  in 
addition,  give  rise  to  intermodulation  effects  producing  sum  and 
difference  frequencies.  In  general,  the  intensity  in  the  region 
in  which  the  radiation  from  different  sirens  can  interact  will 
be  so  much  lower  than  the  intensity  in  the  throat  of  the  siren 
horn  that  these  effects  will  also  be  small  compared  with  harmonic 
distortion  generated  in  the  sirens. 
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PROPAGATION  OF  HIGH  INTENSITY  SOUND  WAVES 


1.  INTRODUCTION 

The  Air  Force  Flight  Dynamics  Laboratory  operates  a  large 
facility  for  testing  airplane  structures  in  intense  sound  fields. 
Very  high  intensity  sound  fields  are  produced  by  banks  of  sirens, 
the  level  reaching  up  to  about  170  db  in  the  vicinity  of  the 
mouth  of  the  siren.  The  question  arose  whether  this  high  inten¬ 
sity  modifies  the  low  amplitude  assumptions  normally  made  when 
computing  sound  pressure  levels  and  whether  the  nonlinear  effects 
of  the  medium,  that  is,  the  air  in  the  facility  will  introduce 
effects  not  accounted  for  by  the  normal  small  amplitude  theory. 

A  project  was  set  up  with  IIT  Research  Institute  under  Contract 
No.  AF  33 (615)- 3436  to  study  possible  effects  due  to  nonlinear¬ 
ities  of  the  air  in  the  facility  as  normally  operated.  Effort 
was  predominantly  directed  toward  computation  of  the  magnitude 
of  nonlinear  effects  from  equations  obtained  from  the  literature. 
In  addition,  where  necessary,  derivations  were  directed  toward 
establishment  of  possible  effects  not  found  in  the  open  litera¬ 
ture.  A  small  program  was  also  run  on  determining  experimentally 
the  effects  of  high  intensity  broadband  noise  as  measured  in  an 
acoustic  tube.  These  measurements  were  further  augmented  by  a 
few  measurements  made  by  the  Air  Force  Flight  Dynamics  Laboratory 
in  their  own  facility  with  two  of  the  sirens  running.  Waveforms 
observed  in  the  AFFDL  facility  are  shown  in  Appendix  A. 
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2.  PHYSICAL  DESCRIPTION  OF  THE  HIGH  INTENSITY  EFFECTS 


For  the  small  amplitude  acoustic  approximation  of  the 
propagation  of  pressure  disturbances  through  a  fluid  it  is 
assumed  that  the  pressure  fluctuation  is  small  compared  to  the 
ambient  pressure  and  that  the  instantaneous  particle  velocity 
is  small  compared  with  the  propagation  velocity  of  the  disturb¬ 
ance.  These  assumptions  relate  to  the  equation  of  propagation: 

c2p"  =  p  (2.1) 

The  second  derivative  of  the  pressure  with  respect  to  distance, 
p",  is  related  to  the  second  derivative  with  respect  to  time,  *p, 
by  the  propagation  velocity  squared.  The  most  predominant  effect 
for  moderate  finite  amplitudes  is  the  influence  of  the  instan¬ 
taneous  particle  velocity  on  the  sound  velocity.  The  sound  velo¬ 
city  is  no  longer  a  constant  but  is  the  sum  of  the  instantaneous 
particle  velocity  and  the  local  sound  velocity.  In  addition, 
the  local  sound  velocity  increases  in  the  regions  of  high  pres¬ 
sure  due  to  the  adiabatic  heating  of  the  medium  and  decreases 
due  to  the  adiabatic  cooling  in  the  low  pressure  region.  In 
this  manner  the  sound  velocity,  c,  is  no  longer  a  constant,  cQ, 
but  is  related  to  instantaneous  particle  velocity,  U,  by  Eq.(2.2). 

c  =  c  +  f3U  (2.2) 

o 

where  p  -  (7  +  l)/2  and  7  is  the  ratio  of  specific  heat  at 

constant  pressure  to  the  specific  heat  at  constant  volume. 

2.1  DISTORTION 

The  higher  propagation  velocity  in  regions  of  positive 
pressures  causes  these  regions  to  overtake  the  regions  of  low 
pressures  causing  a  distortion  or  change  of  waveform  as  the 
pressure  disturbances  are  propagated  through  the  medium.  If  we 
start  with  a  pure  sinusoidal  pressure  fluctuation,  as  shown  in 
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Fig.  2.1(a),  the  overtaking  of  the  high  pressure  regions  will 
cause  a  change  in  the  waveform,  as  shown  in  Fig.  2.1(b).  Solu¬ 
tion  of  the  simple  plane  wave  propagation  equation,  Eq.  (2.1), 
with  c  not  a  constant  leads  to  a  rather  physically  impossible 
situation,  shown  by  the  dotted  line  in  Fig.  2.1(c),  where  mathe¬ 
matically  a  double  valued  pressure  exists  when  the  positive 
pressure  displacement  has  overtaken  the  negative  pressure  dis¬ 
placement.  This  difficulty  arises  only  when  no  loss  mechanism 
exists  in  the  medium.  However,  in  spherical  spreading  or  in  the 
tapered  section  of  the  horn,  this  difficulty  is  not  important 
because  the  drop  of  intensity  with  propagation  eliminates  this 
difficulty.  However,  in  plane  wave  propagation  if  there  is  no 
loss  mechanism  as  postulated,  the  mathematical  solution  is  a 
physical  impossibility  and  a  loss  mechanism  must  be  introduced 
to  the  equations.  Ultimately  the  waveform  assumes  the  shape 
shown  by  the  solid  line  in  Fig.  2.1(c),  in  which  a  shock  wave 
forms,  indicated  by  the  vertical  lines  in  a  pressure- space  or 
pressure- time  diagram.  In  summary  then,  the  sinusoidal  disturb¬ 
ance  of  pressure  propagating  with  large  amplitude  eventually 
changes  its  shape  to  a  sawtooth  form  as  shown  by  the  solid  line 
in  Fig.  2.1(c)  and  remains  fairly  stable  in  this  form,  decaying 
due  to  the  essential  loss  mechanisms  in  the  medium, 

2 . 2  INTERMODULATION 

There  is  an  additional  effect  on  propagation  when  sound 
waves  of  two  different  frequencies  (f^  and  f 2 )  interact  in  a 
region  where  the  medium  is  not  linear  and  where  amplitudes  are 
sufficiently  large  that  nonlinear  effects  must  be  taken  into 
account.  The  high  frequencies,  f^,  move  at  a  higher  velocity 
when  superimposed  with  a  pressure  peak  of  the  low  frequency,  f2» 
than  when  superimposed  on  a  pressure  trough  of  the  low  frequency. 
This  causes  a  frequency  modulation  of  the  high  frequency  sound 
wave  as  shown  in  Fig.  2.2.  Frequency  modulation,  of  course, 
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r 


Fig.  2.1.  Distortion  Generated  at  Large  Amplitudes. 


U 


Fig.  2.2.  The  High  Frequency  Wave  (b)  Superimposed  on  the  Large 
Amplitude  Low  Frequency  Wave  Travels  with  Periodically 
Varying  Velocity  Resulting  in  a  Frequency  Modulation 
as  at  (c). 
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can  be  written  in  terms  of  mean  frequency,  f^,  with  side  bands 
spaced  around  this  mean  frequency  at  intervals  of  the  modulating 
low  frequency  f^  +  nf2*  This  effect  occurs  only  at  points  out¬ 
side  the  siren  in  the  room  itself.  The  most  significant  region 
will  be  the  region  of  spherical  spreading.  However,  for  simpli¬ 
city  a  plane  wave  interaction  series  has  been  developed.  This 
interaction  occurs  at  considerably  lower  intensity  than  the  dis¬ 
tortion  producing  intensity  in  the  throat  of  the  horn.  For  this 
reason  this  effect  is  considerably  less  than  the  distortion 
effect  described  above. 

2.3  BASIS  FOR  THE  MATHEMATICAL  FORMULATION 

Sound  propagation  is,  for  convenience,  subdivided  into 
three  distinct  regions: 

1.  Exponential  spreading  inside  the  horn.  This 
region  of  propagation  is  treated  in  Section  3. 

2.  Spherical  spreading  at  considerable  distances 

from  the  sources.  The  wave  front  will  be  a 

reasonably  good  approximation  to  a  sphere  at 

2  2 

distances  greater  than  about  a  /X ,  where  a 
is  the  area  of  the  horn  mouth  and  X  is  the  wave¬ 
length.  This  region  is  discussed  in  more 
detail  in  Section  4. 

3.  There  is  a  short  intermediate  region  immediately 
outside  the  sirens  where  the  wave  propagates 
essentially  as  a  plane  wave.  This  region  is 
discussed  in  Section  5. 

Interactions  of  waves  of  different  frequencies  are  dis¬ 
cussed  in  Section  6  and  effects  of  reflection  in  Section  7.  The 
effects  of  finite  bandwidth  are  discussed  in  Section  8  and  some 
experimental  results  on  these  are  given  in  Section  9.  This  is 
followed  by  a  brief  discussion  of  the  effects  of  the  wall  treat¬ 
ment  in  the  room  in  Section  10.  Recommendations  and  conclusions 
are  given  in  Sections  11  and  12. 
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To  compute  the  local  sound  velocity  c  in  terms  of  the 
small  amplitude  sound  velocity  cq  and  the  particle  velocity  U 
or  instantaneous  acoustic  pressure  p,  and  ambient  pressure  pQ: 


Hence 


co  =  ^o/po  ’ 

C1  =  +  u 

1/2 


(2.3) 


Pi  £o 
pl  Po 


U, 


Substituting  the  adiabatic  relation 

PlPl'7  *  PoPo'7  * 
fl  .  P1(7'1)/27+^l 


(2.4) 


If  the  pressure  amplitude  is  still  sufficiently  small  that 
the  binomial  expansion  may  be  limited  to  the  first  term 


Pi 


(y-l)/27 


Pi  *  Po 


+  1 


-  1  + 


J 


(7-D/27 


7-1  Pi  -  P( 


(2.5) 


27 


Now  the  excess  pressure,  p  =  Px"P0»  is  related  to  the 
particle  velocity  by  the  medium  impedance  pcQ. 


or 


ULPC0  =  P 


—  -  Vo'cA 

P~ 


7U,c 


1^0 


7'U- 


(2.6) 


(2.7) 
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Substituting  Eq.  (2.7)  in  Eq .  (2.5)  and  this  in  Eq.  (2.4)  gives 
*  co  +  eu  (2.8) 

where  £  =  (-y  +  l)/2. 

This  is  the  nonlinear  effect  valid  for  moderate  sound 
pressures  (p^-pQ<«po)  with  which  all  treatments  of  nonlinear 
acoustics  start. 

We  further  compute  the  conditions  necessary  for  the  form¬ 
ation  of  shock  waves.  Fig.  2.1(c).  To  treat  three  separate 
regions  (plane  waves,  spherical  waves,  and  the  exponential  expan¬ 
sion  within  the  horn),  we  start  off  with  a  general  change  of 
sound  pressure  with  distance. 

p  =  R(x)  pm  sin  co(t  -  x/c) 

It  is  possible  to  substitute  various  specific  functions  for  the 
different  propagation  regions.  For  plane  waves,  neglecting 
absorption: 

R(x)  =  1  .  (2.9) 

If  absorption  is  taken  into  account: 

R(x)  -  e*ax  ;  (2.10) 

in  general,  the  absorption  coefficient  a  is  a  function  of  fre¬ 
quency.  In  the  development  in  Section  5,  the  absorption  computed 

on  the  bases  of  heat  conduction  and  viscosity  is  used.  The  re¬ 
sults  there  presented  are  valid  for  dry  air.  The  presence  of 
water  vapor  introduces  an  additional  variable  leading  to  higher 
absorption  especially  in  the  upper  frequency  region  under  con¬ 
sideration. 

In  the  region  within  an  exponential  horn  having  a  flare 
nix 

constant  m:  Sx  =  S  e  the  sound  pressure  function  becomes: 

R(x)  =  e_mx/2  .  (2.11) 
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In  the  region  of  spherical  spreading  the  intensity  de¬ 
creases  as  the  inverse  square  of  the  distance  and  the  sound 
pressure  decreases  as  the  inverse  first  power; 

R(x)  =  r^/x  ,  (2.12) 

Both  within  the  horn  and  where  sound  radiation  diverges 
spherically,  the  intensity  reduction  due  to  the  increasing  area 
of  the  wave  front  greatly  exceeds  the  internal  losses.  The 
latter  are,  therefore,  of  paramount  importance  only  for  plane 
waves . 

The  condition  for  the  formation  of  shock  waves  is  then 
expressed  mathematically  as  the  rate  of  change  of  the  pressure 
with  distance  becoming  very  large 


dp/dx  oo  .  (2.13-)^ 

We  utilize  the  relations  between  particle  velocity,  sound 
velocity,  and  pressure: 

c  =  cQ  +  3U  ,  (2.8) 

p  =  UpcQ  ,  (2.6) 

dp  =  pcQdU  =  (pcQ/p)  dc  .  (2. 14) 

Furthermore,  the  propagation  of  an  initially  simple  har¬ 
monic  wave  of  pressure  amplitude  pm  is 

p  =  pm  R(x)  sin  <d  (t-x/c)  .  (2.15) 

dx + 

dx  +  dp  ,  (2.16) 

\5c/x  Pco 
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When 


dx 


1 


(2.17) 


dp/dx  jo 

and  shock  waves  form. 

Now  differentiate  Eq .  (2.15): 


(2.18) 


— £  I  =  p  R(x)  cos  [a)  (t-x/c)](xu)/c2) 

.*4 


(2.19) 


The  largest  possible  value  of  a  cosine  function  is  unity, 
so  that  by  substituting  Eqs.  (2.18)  in  (2.19)  we  conclude  that 
shock  conditions  begin  to  appear  when 


Pm  R(*) 


XO) 

C  PCo 


>  1 


(2.20) 


The  sound  velocity  is  assigned  the  mean  value: 


c-c  .  (2.21) 

O 

For  the  plane  wave  condition  without  attenuation  by 
Eq.  (2.9):  R(x)  =  1,  so  that  substituting  Eq.  (2.9)  and  (2.21) 
in  Eq .  (2.20): 

nr  yxy  C  o  P  1 

x  >  -2-  =  °  0  =  6.3  x  103  —  —  (cm]  (2.22) 

C  '  Pm  Pt„2,,f|3  Pm  f 

This  is  the  shock  formation  distance  . 

Table  2.1  gives  the  shock  formation  distance  for  plane 
wave  propagation  in  terms  of  frequency  in  kHz  and  sound  pressure 
level  20  log vT72  pm/0.0002.  The  plane  wave  region  a2/\  is  given 
for  comparison,  for  the  particular  sirens  a  =  1  ft. 
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TABLE  2.1 


SHOCK  FORMATION  DISTANCE  IN  FEET 


Sound  Pressure 
Level,  db 

Frequency, 

kHz 

Pm  0" 

dynes/ cm^  0.1  0.3  1 

3 

10 

171 

I05 

21 

7 

2.1 

0.7 

0.21 

165 

5  x  104 

42 

14 

4.2 

1.4 

0.42 

15  7 

2  x  104 

105 

35 

10.5 

3.5 

1.05 

151 

IQ4 

210 

70 

21 

7 

2.1 

(a  =  1  ft) 

a2/X 

0.1 

0.3 

1 

3.3 

10 

From  this  table  it  appears  that  below  1  kHz  the  plane 
wave  generated  at  the  mouth  of  a  1  ft  square  source  would  start 
to  spread  and  become  a  spherical  wave  before  shock  waves  form. 

On  the  other  hand,  if  horns  are  redesigned  so  that  they  can  pro¬ 
duce  reasonable  pure  wave  forms  at  over  3  kHz,  shocks  may  form 
in  front  of  the  siren. 

For  the  condition  of  propagation  in  a  horn,  Eqs.  (2.11) 
and  (2. 18)  give: 


P^oace-1"1^2  (3 / pc  >  1 


m 


o  - 


/  9 

Now  the  quantity  xe  '  has  a  maximum  value  2 /me  when 

x  =  2/m.  This  means  either  shocks  form  for  x  /  2/m  or  they  will 

never  form.  (Note:  at  the  point  x  =  2/m  the  cross-sectional 

_ 2 

area  of  the  horn  is  e  (=  7.36)  times  the  area  of  the  throat.) 
The  maximum  peak  pressure  in  the  throat  must  be 
3 


Pc, 


me 


P  > 
rm  - 


(2.23) 
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Using  the  concept  of  the  cut  off  frequency,  f  ,  below 
which  the  horn  is  a  very  ineffective  radiator 


gives 


f  =  me  /4tt 
c  o 


pc  f  7P  e  f 

p  ^  ^  e  =  — _£ 

P  f  p  f 

=  —  x  3.16  x  10^  dynes  cm  ^ 


(2.24) 


(2.25) 


Table  2.2  gives  the  sound  pressure  level  at  the  throat  of 
the  horn  for  shock  wave  production  as  a  function  of  the  ratio  of 
frequency,  f,  to  the  horn  cut  off  frequency,  f  . 


TABLE  2.2 

SOUND  PRESSURE  LEVEL  AT  THE  THROAT  OF  THE  HORN 


f/fc  1  2  5  10  20  50 

20  log[^/"0.5  pm/0, 0002]  201  198  194  191  188  184  db 


R(x)  =  rQ/x  and 


For  spherical  radiation,  Eq .  (2.12) 

Eq .  (2.20)  yield  for  the  condition  for  shock  formation 

3 


Pm  ^ 
m  — 


Pc, 


Poor 


7P  c 
,ro  o 

pair 


(2.26) 


6.2  in5  .  -2 
-  x  10  dynes  cm 

Fr 


(2.27) 


where  F  is  the  frequency  in  kHz  and  rQ  the  radius  of  the  wave 
front  where  the  sound  pressure  is  p  . 

It  is  now  necessary  to  relate  the  sound  pressure  level  at 
some  distance  ro  where  propagation  has  developed  to  a  spherical 
form  to  the  sound  pressure,  p  ,  at  the  mouth  of  the  horn. 

O 
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If  losses  are  neglected,  then  the  total  sound  power  through 
a  hemisphere  of  radius  rQ  is  equal  to  the  sound  power  from  the 
mouth  of  the  horn. 


m 


2irr 


2p  2 

o 


1 


(2.28) 


g 


p  r  J  2  rr/S 
Mn  o  T  m 


So  that  conditions  for  shock  waves  forming  in  the  region 
of  spherical  spreading  becomes 


>P0C0  -T** 

Pg  = 

2 

Substituting  for  the  area  of  the  mouth  (1  ft  ): 


(2.29) 


S  =  (30. 4)2  cm2  (2.30) 

-  m 

p  >  5.45  x  10^ /F  dynes  cm  2 

g  " 


Table  2.3  gives  the  sound  pressure  level  L  in  db  re 
^  ?  § 

0.0002  dynes  cm'  at  the  mouth  of  a  horn  as  a  function  of  fre^ 

quency  in  kHz  below  which  no  shock  formation  in  the  region  of 

spherical  spreading  is  expected. 


TABLE  2.3 

SOUND  PRESSURE  LEVELS  AT  THE  MOUTH  OF  THE  HORN  ABOVE  WHICH 
SHOCK  WAVES  WILL  FORM  IN  THE  REGION  OF  SPHERICAL  SPREADING 


F,  kHz  0.1  0.2  0.5  1  2  5  10 

L  ,  db  196  193  189  186  183  179  176 

g 
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Comparison  of  this  table  with  Table  2.2  shows  that  shock 
formation  in  the  region  of  spherical  spreading  will  not  occur 
with  sirens  of  presently  available  power.  Thus,  if  sirens  are 
modified  to  have  a  higher  cut  off  frequency  so  that  the  forma¬ 
tion  of  shock  waves  in  the  horns  is  avoided,  then  shock  forma¬ 
tion  in  the  short  plane  wave  region  immediately  in  front  of  the 
sirens  is  the  most  important  consideration.  Table  2.3  indicates 
that  shock  formation  in  the  region  of  spherical  spreading  need 
not  be  considered.  A  more  precise  treatment  of  shock  formation 
under  conditions  of  spherical  spreading  is  considered  in  Sec¬ 
tion  4.  The  assumption  that  radiation  immediately  in  front  of 
the  horn  approaches  plane  waves  is  a  reasonably  good  assumption 
for  the  horns  of  current  design,  having  a  large  center  core. 

If  horns  with  higher  cut  off  frequencies  are  designed,  this 
assumption  should  be  reviewed. 

2.4  GROWTH  OF  DISTORTION  IN  PLANE  WAVE  PROPAGATION 
OF  INITIALLY  SIMPLE  HARMONIC  WAVES 

The  growth  of  harmonics  in  plane  waves  has  been  computed 

( 1 )  (2) 

taking  full  account  of  absorption  by  Blackstock  and  by  Cook, 

The  method  of  computation  was  different  and  Blackstock  states 
that  his  results  were  sufficiently  close  to  those  of  Cook  that 
he  found  it  not  worth  while  reporting  the  results.  In  both 
cases  digital  computer  techniques  were  used.  In  both  cases  the 
rate  of  growth  of  harmonics  was  computed  as  a  function  of  the 
product  of  absorption  coefficient  of  the  fundamental  (Blackstock: 
a;  Cook:  a  )  and  the  shock  formation  distance  (Cook:  L  *  x 
in  Section  2.4,  Eq .  2.22  of  this  report).  Blackstock  calls  the 
function  aQL  =  T~ * .  Blackstock  extended  computation  over  the 
range  0.2<r<100.  Cook  showed  that,  prior  to  the  formation  of 
a  shock  the  growth  of  harmonics  for  aoL<^0.01,  (T)>100),  (Ref.  2, 
Fig.  9.11),  the  function  of  the  rate  of  growth  of  harmonics  is 
substantially  the  same  as  that  computed  much  more  simply  indepen¬ 
dently  by  Fubini-Ghiron  by  Keck  and  Beyer^^,  and  by 

Hargrove ^ ^  for  the  condition  aQ  =  0. 
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Inserting  into  the  expression  for  ax  values  to  be  ex- 
pected  in  the  operating  range  of  the  facility  under  consideration 


a 


.2/  3 

VO)  /c 


(2.31) 


where 


Hence , 


v'  =  v  [1  +  T)'/2ri  +  (7-l)/2Pr] 
v  -  kinematic  viscosity  —  0.15  for  air  at  70°F 
T)'/t)  =  ratio  of  dilations  to  shear  viscosity 
=  0.8  for  air 

P  *  Prandtl  number  =  0.71 

r 

using  Eq .  (2.22) 


ax 


c 


7P  2?rf 

I  *0 

p  £c  2 
Km  o 


(2.32) 


0.251  x  1.4  x  1.01  x  106  x  6.28  f_ 

1.2  x  (3.44)*  x  10s  p^ 

»  1.58  x  10"3  f/pm  (2.33) 

Now  in  the  problem  under  consideration  f  10^  and  p  ^  10^ 

(151  db)  (for  lower  sound  pressures  the  nonlinear  effects  are 
small);  so  that  in  the  range  of  interest  for  the  facility 
axc  < 0. 01 « 1.  Therefore,  in  the  region  prior  to  the  formation 
of  a  shock  the  attenuation  is  negligible  and  the  simple  expres¬ 
sion  for  the  intensity  of  harmonics  of  Fubini-Ghiron  ,  Keck 
and  Beyer ,  or  Hargrove ^  is  accurate  to  better  than  about 
1  db. 


p 

03 

“  Pm  12  bn  sin  a)<t-x/co> 

(2.34a) 

n=l 

bn 

=  (-1)1*1  J  (n6 ) 

n6 

(2.34b) 

6 

"  x/xc 

(2.34c) 

J  (  )  is  the  Bessel  function  of  order  n. 
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Beyond  the  shock  formation  distance,  Blacks  tock^ 
gives  the  attenuation,  ag ,  of  the  fundamental  as 

a  =  (d  +  1 ) / 2  (2.35) 

s 

The  amplitudes  of  the  harmonics  in  this  region  have 
not  been  computed  in  detail.  However,  for  distances  consider¬ 
ably  in  excess  of  the  shock  formation  distance,  the  wave  may 
be  assumed  to  be  a  sawtooth. 


■pHe  Fourier  series  for  a  sawtooth  wave  is; 


oo 


(2.36a) 

(2.36b) 


Combining  these  last  three  equations,  Eq.  (2.34)  for 
x/xc<l,  Eq.  (2.35)  for  x/xc>2,  and  Eq.  (2.36)  for  x/xc »  1 
(x/x  ^  5  say) ,  a  plot  of  the  amplitudes  of  the  first  four  har 
monies  is  obtained  and  shown  in  Fig .  2.3. 

The  solid  lines  show  the  region  in  which  formulae  are 
valid,  the  dotted  lines  indicate  estimated  interpolations. 

Some  caution  is  needed  in  the  region  x/xc^>)>  1.  Formula 
(2.35)  applies  only  in  the  region  ax  =  0.  For  finite  attenua¬ 
tion,  the  results  of  Blackstock's  computer  program  are  needed. 
Values  are  indicated  by  dashed  lines  in  Fig,  2.3  for  various 
values  of  axc  for  the  fundamental  only.  Again,  from  Eq ,  (2,33) 
ax  is  quite  small  for  all  cases  of  interest.  The  solid  lines 
in  Fig.  2.3  are,  therefore,  reliable  within  the  limits  of  the 
most  sophisticated  computations  available. 

It  should  be  remembered  that  the  theoretical  treatments 
consider  only  the  classical  attenuation.  The  presence  of  water 
vapor  in  the  air  introduces  additional  absorption.  This  will 
attenuate  the  distortion  products  more  than  predicted  by  the 
classical  derivation.  This  will  make  ctx^  somewhat  larger  and 
a  more  complicated  function  of  frequency. 
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Fig.  2.3.  Growth  of  Harmonics  and  Attenuation  of  Plane 
Waves  Due  to  Finite  Amplitude  Effects. 
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2.5  NOTE  ON  NOISE 


All  preceding  sections  have  dealt  with  initially  simple 
harmonic  pressure  fluctuations.  It  is  of  interest  also  to  deter¬ 
mine  the  effects  of  nonlinearities  on  nonperiodic  functions  pro¬ 
duced  by  amplitude  and  frequency  modulation  of  a  siren  output. 

The  taped  program  normally  employed  to  modulate  the  fre¬ 
quency  and  amplitude  of  the  sirens  is  designed  to  simulate  the 
properties  of  band  limited  white  noise. 

The  most  significant  difference  between  band  limited 
noise  and  simple  harmonic  motion  is  the  probability  distribution 
of  the  instantaneous  amplitude.  For  true  white  noise  or  filtered 
white  noise  this  is  a  Gaussian  function  having  a  maximum  at  zero 
and  approaching  zero  for  very  large  amplitudes .  The  exact  Gaus¬ 
sian  function  cannot  be  achieved  either  by  amplification  of  ran¬ 
dom  generators  or  by  modulating  sirens.  In  both  cases  there 
will  always  be  some  finite  amplitude  above  which  the  probability 
amplitude  will  drop  to  zero  much  more  rapidly  than  the  error 
function.  In  the -case  of  modulated  sirens,  this  will  be  the 
maximum  amplitude  the  sirens  are  capable  of  generating  with  the 
modulating  ports  wide  open. 

Similar  limitations  are  in  practice  imposed  on  band 
limited  noise  by  maximum  power  handling  capabilities  of  the 
amplifiers.  The  waveforms  from  the  sirens  may,  therefore,  be 
essentially  similar  to  the  kind  of  waveforms  generated  by  elec¬ 
tronic  means  differing  from  theoretical  Gaussian  noise  in  a 
minor  manner  in  the  region  of  rarely  occurring  very  large  ampli¬ 
tudes.  All  the  probability  functions  have  a  maximum  value  for 
zero  amplitude  decaying  along  the  familiar  bell-shaped  curve  for 
large  amplitudes. 

For  sine  waves  the  amplitude  probability  function  is 
exactly  zero  beyond  the  peak  amplitude;  moreover  this  probability 
amplitude  increases  from  a  minimum  at  zero  amplitude  to  a 
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cusp- shaped  highest  value  at  the  peak  amplitude.  Since  non¬ 
linear  effects  will  depend  on  the  amount  of  time  spent  at  large 
amplitudes  and  the  values  of  this  largest  amplitude,  there  are 
two  opposing  factors.  For  the  same  root  mean  square  amplitude 
the  noise  function  will  be  of  small  amplitude  a  much  larger 
percentage  of  the  time  than  the  simple  harmonic  function.  On 
the  other  hand  considerably  higher  amplitudes  are  reached. 

It  would  be  an  interesting  exercise  to  attempt  a  compu¬ 
tation  on  the  effect  of  this  different  amplitude  distribution 
for  the  two  types  of  functions  considered.  This  has  not  been 
done  within  the  time  and  budget  limitations  of  the  current  pro¬ 
gram.  Instead,  an  experimental  program  was  set  up  to  measure 
the  distortion  generated  by  plane-wave  band- limited  noise  sig¬ 
nals.  The  experimental  program  is  described  in  Section  9. 

The  growth  of  the  measured  "harmonics"  resembles  the 
growth  of  harmonics  of  pure  tones  of  a  different  root  mean 
square  amplitude  as  calculated  from  Eq .  (2.34).  By  means  of 
a  curve  fitting  technique,  it  is  possible  to  estimate  the  dif¬ 
ference  between  the  root  mean  square  pressure  of  noise  and 
simple  harmonic  functions  which  give  rise  to  the  same  amount 
of  signal  distortion. 

The  data  indicate  that  in  the  region  tested  the  occa¬ 
sional  very  large  amplitudes  of  the  noise  signal  outweigh  their 
rarer  occurrence  and  the  same  distortion  is  produced  at  a  lower 
RMS  level  of  noise  than  required  with  pure  tones.  However,  the 
data  are  insufficient  to  predict  this  difference  for  other  fre¬ 
quencies,  RMS  sound  pressures,  or  band  widths.  Further  work 
along  this  line  would  evidently  be  desirable. 
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3. 


DISTORTION  IN  THE  HORN 


3.1  MODERATE  AMPLITUDE- -LOW  FREQUENCY 


The  propagation  inside  the  horn  of  the  siren  occurs  at 
much  higher  intensities  than  elsewhere  and  nonlinear  effects 
manifest  themselves  much  more  strongly  in  this  region. 

The  computations  quoted  by  Olson^^  carried  out  by 
Thuras,  Jenkins  and  O'Neil^  show  that  for  moderate  amplitudes 
the  relative  amplitude  of  the  second  harmonic  p2  is  given  in 
terms  of  the  sound  pressure  at  the  throat  p^  by: 


,—rr  .  ,  ,  -mx/2 

h  =  ii_  Hi  h.  '  _  (3.D 

?1  2  7  PQ  tn  c 

where  m  is  the  flare  constant  of  an  exponential  horn.  We  sub¬ 
stitute  in  this  equation  the  relative  areas  of  the  mouth,  Sm, 
and  the  throat,  St,  of  the  horn,  and  replace  the  flare  constant, 
m,  by  the  horn  cut  off  frequency,  fc  =  mc/4rr.  Further,  let  us 
express  the  relative  level  of  the  second  harmonic  L^  ~  L^  - 
20  log  p2/p^  at  the  mouth  of  the  horn  in  terms  of  the  sound 
pressure  level,  L^,  re  0.0002  p,bar  at  the  throat 


Po 

7 

ui/mc 


20  log  pt  +  74 
6  "2 

10°  dynes  cm  (1  atmosphere) 
1.4  ;  (y+D/y  =  1-72 

u>/2„,c  =  f/2fc 


=  L„  +  20  log  x  1.72  x  106  +  log  (1  -J  5,/Sm)  +  log  — 

4  fc 

-  Lt  +  log  f/fc  +  log  (1  -yj  St/SB)  -  198  (3.2) 

This  function  is  plotted  for  Sfc Sm  in  Fig.  3.1. 
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Fig.  3.1.  Level  of  Second  Harmonic  (L2)  Relative  to  the  Fundamental 
(Li)  as  a  Function  of  the  Sound  Pressure  Level  (L^ )  in 
the  Throat  of  the  Horn.  The  parameter  is  the  ratio  of 
the  frequency,  f,  of  the  sound,  to  the  horn  cutoff  fre- 


The  derivation  assumes  that  the  particle  velocity  is 
small  (though  not  as  small  as  for  the  very  small  linear  theory). 
Experimental  confirmation  of  the  theory  is  presented  for  moder¬ 
ate  amplitudes. The  validity  may  be  expected  to  be  reasonably 
good  up  to  L2  -  L-^  —  -15  db, 

3.2  LARGER  AMPLITUDES- -SHOCK  WAVES 

As  the  distortion  increases,  a  shock  wave  forms  eventually 
and  once  a  sawtooth  wave  form  has  been  established,  this  wave 
form  retains  a  reasonably  constant  shape.  The  Fourier  transform 
of  a  sawtooth  wave  form  yields  a  series  of  sine  functions 

OO 

(-l)n+1  ^  1  sin  nm;  (3.3) 

n=l 

The  amplitude  of  the  second  harmonic  will  be  one  half  the 
amplitude  of  the  fundamental.  On  a  sound  pressure  level  scale, 
the  second  harmonic  will  be  6  db  below  the  fundamental.  The 
levels  of  the  first  ten  harmonics  with  respect  to  the  fundamental 
are  given  in  Table  3.1. 


TABLE  3.1 

LEVEL  OF  HARMONICS  IN  A  SAWTOOTH  WAVE 


number  of  harmonic  123  4  5  6  7  8  910 

relative  level,  db  0  -6  -9.5  -12  -14  -16  -17  -18  -19  -20 


The  interpolation  between  the  level  of  the  second  harmonic 
for  a  fully  developed  sawtooth  wave  form  and  the  moderate  ampli¬ 
tude  for  which  theoretical  derivations  have  been  verified  is 
shown  as  a  dotted  line  in  Fig.  3.1.  The  solid  lines  indicate 
the  levels  of  second  harmonic  distortion  computed  to  the  extent 
to  which  the  theory  is  valid.  The  parameter  on  each  line  gives  the 


22 


value  of  the  fundamental  frequency  with  respect  to  the  cut  off 
frequency.  The  maximum  value  of  the  peak  sound  pressure  in  the 
throat  of  the  horn  for  pure  sinusoidal  excitation  is  one  atmos¬ 
phere  (assuming  a  one  atmosphere  mean  pressure).  The  maximum 
RMS  pressure  is,  therefore,  l/VT  atmosphere.  Maximum  possible 
sinusoidal  sound  pressure  level  is,  therefore,  20  log  (0.713  x  10 
0.0002)  =191  db.  Neglecting  losses,  the  intensity  at  the  mouth 
is  the  product  of  the  intensity  at  the  throat  and  the  area  ratio: 


6 


/ 


Expressed  in  sound  pressure  level  notation 


L  -  L  =10  log  S  /Sf. 
ms  °  m  t 


(3.4) 


The  throat  levels  shown  on  the  abscissa  of  Fig.  3.1  may 
be  converted  to  sound  pressure  levels  measured  at  the  mouth  by 
subtracting  10  log  (Sm/St).  This  formula  is  valid  in  the  region 
where  distortion  is  not  very  severe  (second  harmonic  below  about 
-15  db).  Strong  distortion  products  transfer  an  appreciable 
amount  of  energy  from  the  fundamental  to  the  harmonics  so  that 
in  the  regions  shown  by  dotted  lines  on  the  curve  the  correspond¬ 
ing  sound  pressure  level  at  the  mouth  will  be  somewhat  lower 
than  computed  from  the  area  ratio  relation,  Eq .  (3.4). 

The  sound  pressure  levels  at  which  shock  waves  form  was 
given  in  Table  2.1.  It  may  be  noted  that  these  levels  intersect 
the  extrapolated  low  level  distortion  curves  in  Fig.  3. 1  at  the 
level  L2  -  =  -6  db.  Thus,  there  is  an  additional  point  in 

the  extrapolation  region  which  lends  confidence  to  the  validity 
of  the  extrapolation. 

Further  study  of  Fig.  3.1  reveals  tnat  if  distortion  is 
to  be  minimized,  horns  must  be  used  near  their  cut  off  frequency. 
Experimentally  it  has  been  determined  that  the  distortion  in  the 
present  sirens  is  negligible  at  300  Hz  at  full  output,  but  quite 
severe  with  pronounced  shock  formation  at  1,000  Hz. 
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3.3  UTILIZATION  OF  FACILITY 


The  nonlinear  effects  of  the  propagation  of  high  inten¬ 
sity  sound  waves  may  be  either  reduced  by  siren  redesign  or 
taken  into  account  by  computation  and  calibration  over  a  wide 
range  of  sound  pressure  levels. 

3.3.1  Facility  Modifications 

A  redesign  of  some  of  the  horns  and  operation  of  each 
siren  within  a  more  limited  frequency  range  will  minimize  dis¬ 
tortion.  Thus,  low  frequency  horns  would  be  used  for  generating 
low  frequencies  only.  High  frequency  sounds  would  be  generated 
only  by  the  horns  specifically  designed  for  that  region.  The 
present  location  of  the  stator  motors  limits  the  design  flexi¬ 
bility  in  this  direction  within  budget  limitations. 

At  any  rate,  caution  must  be  exercised  in  replacing  pres¬ 
ent  siren  horns  with  horns  having  a  lower  cut  off  in  order  to 
generate  intense  sounds  at  lower  frequencies.  A  50  Hz  horn  oper¬ 
ating  at  250  Hz  will  generate  as  much  distortion  as  a  200  Hz  horn 
operating  at  1  kHz,  In  other  words,  if  some  horns  are  to  be  re¬ 
placed  with  horns  having  an  extended  low  frequency  range,  this 
replacement  should  be  limited  to  the  minimum  number  necessary. 

3.3.2  Compensation  for  Nonlinear  Effects 

The  present  system  may  be  accepted  as  it  is  with  the 
generation  of  harmonics  at  moderate  and  high  frequencies.  Some 
of  the  sound  power  generated  will  be  diverted  to  harmonics.  If 
this  is  taken  into  account,  there  will  be  no  great  loss  except 
that  at  the  highest  power  levels,  performance  cannot  be  extrapo¬ 
lated  from  measurements  at  low  levels. 
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4.  PROPAGATION  OF  SPHERICAL  WAVES 


In  the  infinitesimal  sound  theory  a  useful  approximation 
of  a  real  sound  source  is  frequently  the  spherical  source.  This 
approximation  is  suitable  when  the  physical  dimensions  of  the 
source  are  small  relative  to  the  wavelength  of  sound  and/or  dis¬ 
tances  from  the  source  are  moderately  large  relative  to  wave¬ 
length.  The  region  in  which  spherical  spreading  is  a  good  approx¬ 
imation  was  given  in  Section  2,  last  line  of  Table  2.1.  If  many 
adjacent  sirens  operate  synchronously,  the  plane  wave  region  ex¬ 
tends  proportionately  further  into  the  room.  Even  then  spherical 
spreading  is  the  predominant  propagation  mode  over  most  of  the 
room. 

Analytic  models  approximating  sources  of  this  type  have 
not  been  incorporated  into  the  problem  of  propagation  of  finite 
amplitude  waves.  To  do  so  requires  the  statement  of  the  funda¬ 
mental  differential  equations  in  an  appropriate  coordinate  sys¬ 
tem  and  the  development  of  an  analytic  statement  of  the  boundary 
conditions  appropriate  to  the  sound  source  and  the  boundaries 
of  its  environment.  This  has  been  done  for  a  piston  radiating 
plane  progressive  waves,  and,  also,  to  a  limited  extent,  for 
spherical  sources.  These  two  cases  represent,  in  a  sense,  limit¬ 
ing  cases  between  which  actual  sources  like  the  sirens  will  be 
expected  to  fall.  The  plane  wave  case  represents  that  extreme 
wherein  loss  of  intensity  in  the  wave  occurs  only  because  of 
loss  to  the  medium. 

4.1  APPLICATION  OF  THE  COMPUTATIONS  BY  LAIRD  AND  ACKERMAN 
TO  THE  SONIC  FATIGUE  FACILITY 

One  approach  to  the  analytic  treatment  of  the  propaga¬ 
tion  of  spherical  waves  of  finite  amplitude  is  the  paper  by  Laird 
and  Ackerman. ^  They  begin  with  the  usual  equations  of  con¬ 
tinuity,  motion,  and  the  adiabatic  equation  of  state  for  a  per¬ 
fect  gas  in  the  Eulerian  formulation.  With  the  additional  assump¬ 
tion  of  spherical  symmetry  about  the  origin,  these  equations  can 
be  stated  as  : 
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0 


(4.1) 


+  U  +  p  —  +  — ^  - 

5t  5r  dr  r 

p  —  +  pU  —  +  ^  -  0  (4.2) 

dt  dr  dr 

-e-  -  M  (4.3) 

P0  W 

where  p  =  density,  p  =  pressure,  and  U  =  particle  velocity. 

The  assumption  that  the  gas  is  non-viscous,  non-heat  con¬ 
ducting,  and  non- heat  radiating  is  a  common  one  in  the  literature; 
and  it  is  generally  assumed  that  the  approximation  is  satisfactory 
during  that  regime  in  wave  propagation  before  shock  waves  begin 
to  form. 


After  several  transformations  of  coordinates,  including 
a  normalization,  Eqs.  (4.1)  and  (4.2)  become 

*S  +  _2_  v  dS  +  (I+S)  +  2v(l+S2  =  0  (4.4) 

7-1  d0  7-I  dx  dx  x 


—  +  v  —  +  — (1+S)  —  =  0  (4.5) 

d0  dw  7- 1  dw 


where  v 
S 


w 

0 

f 


c 

c 


o 


U/co 

<c-co>/co 
(27 r/XQ)  r 

27Tf  t 

frequency 

wavelength  of  sound  in  the  limit  of  infinitesimal 
amplitude 

yj  7p/p  *  velocity  of  sound 

=  undisturbed  speed  of  sound 
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Boundary  conditions  are  specified  in  terms  of  particle  velocity 

at  the  surface  of  a  pulsating  spherical  source  with  equilibrium 

radius  r  ,  in  the  form 
o 

11^  =  U  (rQ  +  £p)  =  Um  exp  j  (2irft-^)  (4.6) 

£  = - —  U  exp  j  (2-rrf  t-  )  (4.7) 

p  27rf  m 

where  Um  is  velocity  amplitude  of  the  source  and  $  is  an  arbi¬ 
trary  phase  angle.  We  note  that  the  boundary  conditions  do  not 
give  the  particle  velocity  at  a  fixed  radius  rQ  but  at  a  vari¬ 
able  radius  r  +  £  where  £  is  the  displacement  of  the  surface 
o  "p  P 

of  the  sphere  from  equilibrium  at  time,  t. 

In  normalized  coordinates,  Eqs.  (4.6)  and  (4.7)  are 


v  =  v  (w  +  K  ,y)  =  v  exp  j  (©-(£)  (4.8) 

p  o  p  m 

Kp  =  -jvm  exp  j  (€i-«()  (4.9) 


where  v  =  h,/c0-  Since  this  formulation  gives  particle  velocity 
at  time  t  at  a  variable  radius  rQ  +  t ^ ,  it  is  difficult  to  apply 
the  boundary  conditions  in  order  to  obtain  a  specific  solution. 
This  is  overcome  by  expanding  particle  velocity,  v  „ ,  in  a  Taylor's 
series  in  terms  of  particle  velocity  at  radius  rQ.  That  is,  we 
obtain  an  approximation  to  what  the  particle  velocity  as  a  func¬ 
tion  of  time  must  be  at  rQ  to  produce  the  given  particle  velocity 

at  r  +  £  at  the  same  time . 

o  p 


Laird  and  Ackerman  solve  the  equations  by  a  perturbation 
method  out  to  the  second  order  approximation.  The  limiting 
assumptions  made  in  this  procedure  include  the  condition  that 
the  source  circumference  shall  not  be  significantly  smaller  than 
one  wavelength,  i.e.  ,  wQ  )>)>  1 .  Since  the  horn  opening  to  the 
sirens  is  of  the  order  of  1  ft  in  radius,  it  seems  reasonable  to 
assume  that  this  condition  will  be  met  down  to  frequencies  of 
200  Hz  (wavelength  approximately  5.65  ft).  The  total  solution 


out  to  second  order  has  the  form 
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v(w,0)  =  A^^exp  j(O-w)  4-  +  ^^exp  j2(0-w)  (4.10) 


S(w,9)  =  [Afi^^exp  j(0-w)  +  CSq^2^  +  <52^2^exp  j2  (0- w)  ]  (4. 11) 

2 


w 


where  A  =  v_ 


m  (vAi)1'2 


vm  =  normalized  particle  velocity  amplitude  at  source 
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The  actual  solutions  are  to  be  the  real  part  of  these  complex 

7  8) 

quantities.  A  computer  program  was  used  to  obtain  numerica 
results  in  a  reasonable  time. 

Expanding  the  perfect  gas  law  into  a  Taylor  series  pro¬ 
vides  an  analytic  expression  for  pressure,  p,  in  terms  of  the 
solutions  for  S,  the  normalized  local  velocity  of  sound.  Sub¬ 
stituting  in  the  solution  for  S  up  to  second  order  gives 


(4.12) 


The  computer  calculation  gave  v,  S,  and  p  as  functions  of  nor¬ 
malized  distance,  w  =  (2tt/\  )r  and  v  =  U  /c  .  The  following 
tables  give  the  range  of  these  parameters  of  interest  to  the 
Sonic  Fatigue  Facility  program. 


TABLE  4.1 


NORMALIZING  FACTOR  VERSUS  FREQUENCY 


200 

500 

1000 

2000 


0.908 

0.363 

0.182 

0.091 


1.10 

2.75 

5.50 

11.00 
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TABLE  4.2 


NORMALIZED  VELOCITY  AMPLITUDE  VS.  SPL 


Sound  Pressure  Level 
re  0.0002  pbar 


114 

134 

154 

174 


The  computer  program  used  the  parametric  values  of 
wq  *  0.2066,  wq  =  1.0,  and  wQ  =  5.0;  and  vm  ~  0.05,  vm  -  0.1, 
v  =0.2,  and  v  =0.4.  Of  these  values,  w  =  1.0  and  w  =  5.0 

Trt  m  u  v 

together  with  =  0,05  are  of  interest  for  our  case. 

This  conventional  perturbation  solution  fails  to  converge 
at  large  distances  from  the  source*  Therefore  ,  for  large  values 
of  w ,  a  perturbation  solution  was  developed  based  on  the  charac¬ 
teristic  equations  associated  with  the  basic  system  of  differen 
tial  equations,  Eqs.  (4.1),  (4.2) ,  and  (4.3).  The  procedure  is 
based  on  that  discussed  by  Lin. ^  ^  The  solution  was  carried 
through  the  second  order,  but  only  the  first  order  solution  was 
considered  in  discussing  quantitative  results.  When  the  ampli¬ 
tude  is  low  enough  s+v<^<(1,  this  perturbation  of  characteristic 
equations  solution  can  be  used  to  describe  the  formation  of  the 
shock  fronts,  their  growth  and  decay*  When  the  amplitude  is 
high  near  the  source  but  is  appreciably  diminished  by  spherical 
divergence  before  shocks  are  formed,  it  is  convenient  to  use 
the  conventional  perturbation  solution  carried  to  second  or 
higher  order  near  the  source  out  to  some  intermediate  distance 
where  the  amplitude  is  lower  and  then  to  use  the  first  order 
characteristic  perturbation  solution  for  greater  distances. 
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This  was  done  in  the  computer  program  by  Laird  and  Ackerman.  In 
Figs.  4.1  and  4.2  we  have  replotted  their  results  for  two  cases 
applicable  to  the  Sonic  Fatigue  Facility. 


In  Fig.  4.1,  wQ  =  (2Tr/XQ)ro  was  chosen  to  be  1.0.  This 

corresponds  to  a  frequency  of  180  Hz  if  r  ,  the  radius  of  the 

source,  is  taken  as  1  ft.  vm  =  U  /c  was  selected  as  0.1  about 

m  mo 

174  db  SPL.  The  plot  of  sound  pressure  level  versus  distance 
for  the  fundamental  falls  closely  upon  the  curve  for  pure  spheri¬ 
cal  spreading,  i.e.,  6  db  per  doubled  distance.  The  second  har¬ 
monic  appears  prominently  at  only  1  ft  from  the  source,  but  decays 
rapidly.  The  calculation  of  the  value  of  the  normalized  distance 
at  which  the  shock  wave  begins  to  develop,  wg^,  gives  a  value 
y  10^  or  approximately  1,000  ft.  The  expression  for  wg^  is  based 
on  the  characteristic  perturbation  solution 


sh 


7+1  v  w 
m  o 


<^o  > 


w. 


1/2 


(4.13) 


It  has  been  assumed  that  the  source  is  sufficiently  weak  so  that 
the  first  order  perturbation  in  the  characteristics  may  be  used 
to  describe  the  whole  field  and  wq  is  of  the  order  of  magnitude 
of  1. 


The  calculations  indicate  that  for  frequencies  at  the 
lower  end  of  the  spectrum,  spherical  divergence  would  preclude 
the  formation  of  shock  waves  within  the  Sonic  Fatigue  Facility 
due  to  propagation  in  the  air  of  the  room.  This  is  true  for 
even  the  very  highest  sound  pressure  levels  contemplated.  The 
presence  of  higher  harmonics  in  the  sound  field  at  these  fre¬ 
quencies  would  be  due  to  their  generation  inside  the  sirens  them¬ 
selves.  In  other  words,  since  spherical  divergence  character¬ 
izes  the  sound  field  at  low  frequencies,  we  cannot  expect  the 
generation  and  growth  of  the  higher  harmonics. 
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Fig.  4.1.  Growth  of  Second  Harmonic  in  a  Spherical  Wave  of  Finite  Amplitude. 


180 


o 

o 


o 


O 


o 

m 


o 

<r 


o 

to 


o 

CM 


jeq  2000  *  0  3a  qp  aJtnssajj;  punos 


CM 

m 

*h 

Uh| 


33 


Growth  of  Second  Harmonic  in  a  Spherical  Wave  of  Finite  Amplitude. 


Figure  4.2  displays  the  results  of  the  computer  program 

for  a  choice  of  w„  s  5 , 0  and  v  =0.05.  A  choice  of  w  =5 

o  m  o 

corresponds  to  a  frequency  of  about  1,100  cps  if  rQ  is  chosen 
as  1  ft.  The  second  harmonic  grows  to  a  significant  component 
of  the  signal. 


4.2  APPLICATION  OF  RECENT  RUSSIAN  LITERATURE 
TO  THE  SONIC  FATIGUE  FACILITY 


A  more  recent  treatment  of  the  propagation  of  finite 
amplitude  spherical  waves  is  that  of  Naugol'nykh,  Soluyan,  and 
Khoblov. ^ This  paper  represents  an  attempt  to  include  con¬ 
sideration  of  viscosity  and  thermal  conductivity  of  the  medium. 
To  do  this  the  Navier- Stokes  equation  is  formulated  in  the  form 

p  fay  +  u  su 

at  Sr 

where  b  =  4/3  q  +  q'  +  k 


-  p  -  lE  +  b 

[l  -  2,  « 

hr 

r  dr^  r^ 

(4.14) 


c 


V 


c 


PJ 


represents  the  results  of  combining  the  viscosity  terms  in  the 
Navier- Stokes  equations  and  the  heat  terms  in  the  equation  of 
state  taking  into  account  the  fact  that  the  process  is  non- 
adiabatic.  The  total  system  of  equations  to  be  solved  is 

completed  by  the  equation  of  continuity  and  the  approximate 
equation  of  state. 


il£  + 

U  -^  +  p 

au  +  2£U  =  0 

(4.15) 

dr 

^r  r 

P  = 

,  2 
PQ  +  c 

<p-p0>  +2r  ^2  (p-p=)2 

(4.16) 

o 


A  series  of  transformations  and  approximations  leads  to 
the  final  differential  equation  to  be  solved 
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(4.17) 


su  +  u  _  AU  su  =  6  a2u 

Sr  r  St  Sr^ 

where  A  =  ;  5  =  — - -  ;  t  =  t  -  — 

2  3 

2c/  2p  c  J  c 

o  00  o 


For  the  case  of  air,  values  for  these  constants  are: 

b  2  9.12  x  10'6  Newt°"  sec 

raz 

A  =  9.90  x  10" 6  sec2 

meters^ 


6  =  9.22  x  l(f  14  Newton  sec4 

meter^  kgm 

Equation  (4.17)  is  solved  in  the  case  of  large  values 
of  the  reciprocal  dissipation  factor,  R,  i.e. , 


r  =  — a  »  1  ,  i.e. ,  6=0.  (4.18) 

2o;6 

In  our  case 

r  =  *  2.56  x  10®  x  — 

4-rrf  (9.22)  x  10' 14  f 


where  f  is  frequency.  For  all  frequencies  of  interest 
f  <  10,000  Hz,  R»  1,  and  this  solution  should  apply.  In 

Eq.  (4.17)  we  can,  therefore,  drop  the  terra  on  the  right  side 

of  the  equality  sign.  The  resulting  equation  together  with 
the  boundary  conditions  U  =  UQ  sin  (ojt-krQ)  at  the  surface  of 
the  pulsating  medium  r  *  rQ,  leads  to  the  solution 

cut  -  arcs  in  e  -  AcoU  r  log  —  e  (4.19) 

O  O  6 

ro 

i_  Ur 

where  e  =  ■ 
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This  solution  describes  the  propagation  of  the  wave  prior  to 
the  formation  of  the  shock  wave  where  dissipative  process  plays 
no  substantial  role,  ryy  1.  Except  for  the  second  term, 

GU0o>ro  log  (r/rQ),  this  solution  would  be  a  spherically  diverg¬ 
ing  sinusoidal  wave.  The  term  X  =  Go>U  r  log  (r/r  )  results 
in  a  progressive  steeping  of  the  wave  front.  At  any  instant,  t, 
this  factor  increases  logarithmically  with  r.  If  X  is  less  than 
1  this  distortion  is  modest  for  X<^1,  strong  when  X  is  compar¬ 
able  to  1,  and  at  X)>  1,  arr  becomes  multivalued.  The  condition 
that  X  =  1  can  be  used  as  a  criterion  to  define  the  point,  r, 
where  the  shock  wave  begins .  This  leads  to 


l 

;o 


—  *  exp 

r 


kr 


(Ybl)Uc 

2c_ 


=  exp 


1+1  wovmo 


(4.20) 


which  is  equivalent  to  Eq .  (4.13)  of  the  earlier  theory  under 
the  approximation  that  wq  is  about  1, 


As  an  example:  f  =  3000  Hz  and  v  =  U  /c  =  0.05, 
r  moo 

r  =  0.304  m  =  30.4  cm  =  1  ft,  r,  becomes  2.7  r  ,  i.e.,  2.7 
o  *  1  o  ’ 

ft.  The  plane  wave  shock  formation  distance  at  this  intensity 
(171  db  re  0.0002  dynes  cm  from  Table  2.1  would  be  0.7  ft. 
This  shows  the  extent  to  which  spherical  spreading  by  reducing 
intensity  also  reduces  the  distortion.* 


The  two  theoretical  treatments  give  essentially  the 
same  results  except  that  Laird  expresses  his  results  for 

pressure  in  a  form  that  provides  estimates  of  the  pressure 


*Naugol ' nykh,  et  al,  state  that  their  assumption  of  negligible 
attenuation  (R  »  1)  is  valid  only  for  the  region  in  which  dis¬ 
tances  are  not  very  large  compared  with  rG.  The  formula  is, 
therefore,  not  valid  for  low  frequencies  or  moderate  ampli¬ 
tude  when  shock  formation  distances  would  become  very  large. 
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levels  for  each  harmonic  as  a  function  of  radial  distance  from 
the  source  while  Naugol'nykh^  ^  presents  a  much  simpler  rela¬ 
tionship  from  which  the  wave  form  may  be  obtained  as  a  function 
of  distance.  Laird's  procedure  requires,  however,  a  computer 
program  (already  developed  and  used)  to  obtain  results  for  a 
specific  value  of  the  parameters.  This  is  especially  true  if 
we  wish  information  about  harmonics  higher  than  the  second. 
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5.  FINITE  AMPLITUDE  PLANE  ACOUSTIC  WAVES 


The  most  general  differential  equations  involving  wave 
propagation  in  a  gaseous  medium^  ^  are  quite  complex  mathemat¬ 
ically  and  at  the  present  time  are  unsolvable.  However,  by  mak¬ 
ing  suitable  approximations  these  differential  equations  may  be 
linearized  and  the  resulting  simplified  equations  can  be  solved 
to  yield  solutions  applicable  to  "small"  signal  propagation. 

With  the  advent  of  more  powerful  acoustic  sources,  many 
of  the  approximations  made  in  obtaining  the  "small  signal  equa¬ 
tions"  cannot  be  used  if  a  useful  analysis  of  the  propagation 
is  to  be  made.  A  consequence  of  this  limitation  is  that  the 
differential  equations  of  wave  propagation  become  nonlinear. 

Two  methods  of  approach  have  been  used  in  an  attempt  to  solve 
the  nonlinear  problem.  The  first  and  most  common  method 
is  that  utilizing  perturbation  techniques  to  obtain  solutions. 

A  disadvantage  of  this  method  is  that  it  is  relatively  weak  and 
is  tedious  to  use  in  practice.  The  second  method  is  that  of 
transforming  the  equations  into  a  nonlinear  equation  whose  solu¬ 
tion  is  known.  Blacks tock^5  shows  how  the  nonlinear  plane 

wave  equations  may  be  transformed  into  an  equation  called  Burger's 
equation  whose  solution  for  initial  value  problems  is  known. 

The  dependent  variable  is  particle  velocity.  It  is  this  method 
which  will  be  used  in  this  discussion. 

An  approximate  form  of  the  equations  of  motion  in  a  thermo- 
viscous  medium  is  given  by  Blackstock  as: 

co\  -  PcoUUc,  =  ~  v  [Y  +  (7-D/PgJ  Ut,t,  =  v,Ut,t,  (5.1) 

where  x  is  distance,  t'  -  t  -  x/c  ,  t  is  time,  cq  is  small  sig¬ 
nal  sound  speed  in  the  absence  of  dissipation,  v  iv  kinematic 
viscosity,  Y  =  2  +  q'/q  is  the  viscosity  number,  q  and  q'  are 
the  shear  and  dilational  viscosity  coefficients  respectively, 

7  is  the  ratio  of  specific  heats,  and  Pr  is  the  Prandtl  number. 
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The  quantity  p  is  given  by  1/2 (7+I)  for  perfect  gases.  Divid- 

3 

ing  through  by  cq  and  letting  z  =  t'/A,  one  obtains  the  equation 


u  -  uu  =  cu 

x  2  22 


(5.2) 


where 


A  =  p/c 


C  = 


c  v 
o 


2(3' 


Using  the  transformation 


U  =  2C  [low  Cl 


the  equation  reduces  to 


C  =  cc 

^x  zz 


(5.3) 


(5.4) 


where  it  has  been  assumed  that  C  ^  0.  Making  a  further  change 
of  variables  x  -p>  y  =  xc,  one  obtains  the  equation 


C  -  C 

^zz 

which  has  the  initial- value- problem  solution 

„O0  « 


(5.5) 


(15) 


C(y,z)  = 

where 


-±1 

4t r  •'-0 


^(P> 


(y) 

^ (p)  =  C(0,p) 


1/2 


exp 


-  ^z~P^  dp  for  y>0  (5.6) 


4y 


(5.7) 


In  the  case  at  hand 

C(0,p)  =  exp 


^ f 

2C  J0 


U(0,m)  dm 


(5.8) 


This  derivation  was  included  here  as  the  basis  for 
both  Section  6  and  Section  7. 
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6.  SUPERPOSITION  OF  FINITE  AMPLITUDE  WAVES 

Let  us  now  consider  the  case  where  two  plane  progres¬ 
sive  waves  which  are  sinusoidal  at  the  origin  are  superimposed 
throughout  space.  Thus, 


U(o,m)  =  U^  sin  (o^Am)  +  Ug  sin  (o^Am  +  6)  ;  t  )>  0 
U(o  ,m)  »  0  ;  t<0 


(6.1) 


at  the  initial  point  of  superposition.  Evaluating  the  integral, 
Eq.  (5.3),  we  obtain  the  result  that 


C(0,p)  =  B  exp 
where 


-U 


2Ccd^A 


1  u2 

cos  (cu^Ap)  -  -  cos  (o^Ap  +  6)  (6.2) 


2C(X>2^ 


B  =  exp 


U,  U9  cos  <5 

— - —  +  — - 


2Ccd^A  2Ca)2A 


(6.3) 


Thus,  from  Eq .  (5.6) 


Z(y,z)  = 


4=  f 

V  4  Try  Jq 


exp 


-U. 


201^  A 


cos  (oj-^Ap)  + 
2 


-U2  (2-p) 

cos  (aj?Ap  +  6)  -  - 

4y 


2 00^2  A 


dp  + 


B 


\/  4"rry 


f 


exp 


 .(z~p) 

4y 


-  CO 

2 


_U1  -u? 
cos  (a^Ap)  + 


2Ca)-^A 


2Ca>2A 


cos  (ci^Ap  +  6) 


J 


dp 


(6.4) 


The  second  integral  will  vanish  in  the  steady  state,  z  -)>  =>. 


Making  use  of  the  relations 


(16) 


-z  cos  0 


-  E  <-i>n  s 


cos  n9 


n=0 
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n 


1  if  n  =  0  | 

I  n’ 

2  if  n  >  1  1 


an  integer 


l^(z)  =  Modified  Bessel  Function 


we  obtain 


C(y3z)  = 


B 


4iry 


ft£  <- 

•'O  f\R=0 


U 


i)kekIk  " —  cos  (^AP) 


2Cu>1A 


x  I  ^  (-  1) 
>0 


u. 


Hh 


2COJ2A 


cos  (^u^Ap+ifi )  j  exp 


-  (z-p)' 


4y 


dp 


(6.5) 


C(y,z)  = 


B 


CO  00  i*  GO 

A(k,0)  /  cos  (kcu-jAp)  cos  (i^Ap  +  j|6) 
k=0  £=0  •'O 


k=0  i=Q 


*  exp  ^  ^  dp 


where 


U, 


A(kJ)  =  (-l)k+i  e.e  Xk  ( —  t  I, 

K  K  \  2Co)jA 

Making  use  of  the  following  relations 


U, 


ZCo^A 


(6.6) 


(6.7) 


JC°°  \  2  n 

cos  (ap)  exp  "  dp  «  yTkiry  cos  (az)  exp  (-a2y) 

0  4y 

/oo  2 

sin  (bp)  exp  ~  ^z~ P^  dp  =  V  4fry  sin  (bz)  exp  (-b^y) 

4y 


Lim 

z-^<*> 


Together  with  various  trigonometric  relations,  we  obtain  the 
value  of  the  integral  in  Eq.  (6.6)  for  the  steady  state  (z-^>  °°) . 
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/■“  2 
I  cos  (kci'^Ap)  cos  +  6i)  exp  ^Z~  ^  dp 

_  y/'i'rry  cog  [(^^kcL^)  Az  +  £&]  exp  [- (ia^+kco^)2  A2y]  + 


I  «  Lim 

Z-^oo 


cos  [  (/o^-ko)^)  Az  +  6  fc\  exp  [-  (/u^-ka)^)2  A2y]  (6.8) 


Thus 


00  00 


C(y»z)  = 


B 


Y  Z)  A(k-i>  1 


(6.9) 


'^y  i^o  I-o 

and  the  particle  velocity  becomes 

U(y,z)  =  2C  [log  az  (6.10) 

where  z  =  t'/A  and  y  =  xC 

In  the  expanded  form,  the  equation  is 

/  O0  00 

U  =  2C(S  It, 

\k=0  1=0 

(iiu^+ku;^  )  A  exp  [- (j£«J2+ka^)  2A2y  ]  sin  [  (jg^+ka-^ )Az+S/  ]  + 
(/uj2-ka)L)A  exp  [- (ia)'2-k(o1)2A2y]  sin  [  (i^2“k^1)Az+£6]^||^^ 
y;  (-l)k+^eke^  Ik(U1/2Ca)1A)I^  (U2/2Ccd2A)  |cos  [  (^a^+ka^jAz+S £] 

l=o 

0  O 

•  exp  [- (i^2'hulk^  A  yl  +  cos  ( (i  0^2  -ka)^)Az] 


(_  l)k+^+1€k^  ^(U^Ca^A)  1^  (U2/2Ca)2A) 


exp  [- (ia)2-km1)2A2y] 


(6.11) 


For  the  case  of  a  single  source,  we  let  U2  and  5  approach  zero. 
Noting  that 

0  if  1  +  0 

Lim  10  (U„/2Ca>2A)  = 

u2->  0  l  1  if  £  =  0 
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We  find  that  the  particle  velocity  due  to  a  single  source  can 
be  given  as 


U  =  2C 


^  (-l)k+1  ekIk  (U1/2Cm1A)kco1A  exp  (-k2cu12A2y)  • 

00 

^2  C-DkekIk(U1/2CaJlA)  exp  (-k^Vy)  • 

(6.12) 


_k“0 
sin  (kaj^Az) 


cos  (ka^Az) 


k=0 


(It  is  to  be  noted  that  this  solution  is  identical  to  that  given 
by  Blackstock.)  The  result  in  Eq.  (6.12)  may  be  simplified  by 
noting  that  for  large  arguments 


IR(z)  ~  (eZU  2ttz) 


4k2- 1  +  (4k2- 1) (4k2- 9) 
8z  21  (8z)Z 


(6.13) 


(14 ) 

Blackstock  has  shown  that  we  need  use  only  the  first  term 
of  the  asymptotic  expansion  for  IR(z)  when  y  is  large.  This 
is  so  because  in  the  product  of  the  exponential  term  and  the 
term  in  Ik(z)  in  Eq.  (6.12)  the  exponential  term  dominates  the 
Ik(z)  term  when  y  is  large,  i.e.,  this  product  approaches  zero 
as  y  becomes  large.  The  restriction  that  y  be  large  means  that 
the  solution  obtained  from  Eq .  (6.12)  by  using  the  first  term 

of  the  expansion  for  I^(z)  in  Eq.  (6.13)  will  probably  not  be 
valid  near  the  source . 


Making  the  substitution  for  IR(z)  from  Eq .  (6.13)  into 
Eq .  (6.12),  we  obtain 

^  (- l)k+'L(eR)  (kojjA)  exp  [-k2u>i2A2y] 

U  -  2C  — - -  (6.14) 

y]  (~1)R  tR  exp  [~k2oj^2A2y] 
k=0 
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The  quantity 


=  (~l)k  ek  exp  (-k2to^2A2y)  cos  (krn^z) 


k=0 


=  0, 


r^iAz  2  o  i 

1^— —  >  exP  (-ajj^A  y)| 
(16) 


where  0^  is  a  theta  function. 


Now 


(e4)u(u,q)  _  ,  v>  qn 

-  -  4  sin  2nU 


e4(u,q)  — [  l-q 

We  see,  therefore,  that 

k££).  -p  (-^Vy)] 


<1 


U(x,t)  =  2C 


k  (“T-)’  exp  <-“i2A2y>] 


(uj^Az/2),  exp  (-u^2A2y)  ]  (oj^Az/2)  0SjA 
[64  (tOjAz/2),  exp  (-0J12A2y)]  2 


=  2C 


or  simplifying 


2  2 

Eexp  (-tb,  A  yn) 

- n — * -  sin  (nrnAz) 

n=l  1  "  exP  (" 20)^  A  yn) 


2CujlA  "T  ,  , — T7T 


sin  (ncu^Az) 


n_^  sinh  (a>i  A  yn) 


=  2  Cm 


w 

iaE 


sin  [na>1(t  -  x/cD)] 


n=1  sinh  (m12A2xCn) 


(6.15) 


(6.16) 


(6.17) 
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7.  NORMAL  INCIDENCE  REFLECTION  FROM  A  HARD  WALL 


UNDER  SPECIAL  CONDITIONS  (NORMAL  INCIDENCE) 

An  approximate  solution  to  the  problem  of  reflection 
from  a  wall  may  be  obtained  if  the  following  assumptions  are 
made.  Let  a  finite  amplitude  plane  wave  source  be  located  at 
a  point  x  *  0.  Let  a  perfectly  reflecting  wall  be  located  at 
a  point  x  *  p,  far  enough  from  the  source  so  that  the  attenua¬ 
tion  of  the  primary  wave  makes  possible  the  linear  superposi¬ 
tion  of  two  waves  in  this  region  (as  in  the  small  signal  case). 


Under  the  above  conditions,  the  totally  reflecting  wall 
may  be  replaced  by  an  image  source  satisfying  the  boundary  con¬ 
dition  that  the  sum  of  the  particle  velocities  due  to  the  primary 
source  and  the  image  source  is  zero  at  the  position  of  the  wall 
(x  =  p).  From  Eq.  (6.12)  the  primary  source  may  be  given  as 


Uj  *  2C 


sin 


T:  (-l)k+1  ekIk  (U1/2Cai1A)  exp  [- (ku}jA)2Cx]  (kaijA) 


k=0 

!k(Ui  Ac 

~  ~ 

P 


~\  r  « 


co*/=J  T,  (-vk  <v2c“>iA)  • 

’  k=0 


exp  [-(kcu^A)  Cx]  cos 


kuu-,  Ac  ) 

-V"  (t-xM 


(7.1) 


where 


y  =  xC 


z  =  -2-  (t-x/cQ) 


For  the  image  source  the  particle  velocity  satisfying  the 
boundary  conditions  is 
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“  00 


UR  =  2C 


C-l)k+1  ekIk  (U1/2CtD1A)  exp  [  (kujjA)2 (x- 2p)C] 


k=0 


f  ktu,  Ac  /x-2p 
(kcojA)  sin  l — — - ^ 

exp  [  (koj^A)2(x- 2p)C]  cos 
The  composite  field  is  thus 


-  t 


ktu 


L  (-DkekIu(U1/2QUlA) 
k=0 


lAco  /X~2P 

3  V  c~ 


-  t 


(7.2) 


UTotal  UI  +  UR 


(7.3) 
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8.  AUTOCORRELATION  PROPERTIES  OF  A  BAND  OF  FINITE  AMPLITUDE, 


PLANE  WAVE ,  WHITE  NOISE 

The  autocorrelation  function  of  a  signal  may  be  obtained 


(22  ) 

from  the  following  Fourier  transform,' 


J  0 


where  R(t)  is  the  autocorrelation  function  and  S(f)  is  the 
power  spectral  density.  To  an  approximation  (which  is  valid 
for  the  sound  pressure  levels  expected  in  the  RTD  facility) , 
the  sound  power  in  a  plane  wave  is  proportional  to  the  square 
of  the  sound  pressure,  where  pressure  is,  in  turn,  proportional 
to  the  particle  velocity  U.  Thus, 


00 


(8.2) 


where  z  =  constant  =  pc 

o  o 

U  =  RMS  particle  velocity  spectral  density 
s 

The  RMS  velocity  spectral  density  for  the  initially  finite 
amplitude  plane  wave  is , 


U 


s 


(8.3) 


where  t  is  the  period  1  / f . 
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For  the  sinusoidal  plane  wave  the  correlation  function  is 


Ro^Ct)  =  —  —  I  cu2  cos 


7T 


CO  CO 

CD2  COS  (cut)  /J 

Jo  n®l 


1 

T,  2/2. 


sinh  (cu  A  xCn) 


2  2  2  “ 

ZC  A  GD^  COS  ^  ^ 

s  inh2  (cu-i  2A2xCn) 
n-l  i 


6  (oj"ajj)dto 


(8.4) 


TT 


where  6(cu)  is  the  Dirac  Delta  Function.  The  correlation  coef¬ 
ficient  may  be  determined  to  be 


R  <t)  =  R  (r)/ R  (0)  =  cos  cu  r  (8.5) 

“CU^  CO-^  £U  ^  I 

For  the  case  of  a  finite  band  of  noise  cu^  <  cu  <(  cu2  >  Eq. 
(8.2)  becomes 


R(t) 


H 


zC2A2 


-  00 

f  2  ?22 

2C  to  A  cos  a rr 

£ 

i 

sinh2  (cu2A2xCn) 

i  O)-^ 

n=l 

dcu 


7T 


£ 

n=l  cu^ 


CU  COS  CUT 

2,  2.2  ' 


dcu 


(8.6) 


sinh  (cu  A  xCn) 

This  integral  can  be  evaluated  in  the  following  manner.  Let 

(8.7) 


I  = 


f  m2  cos  m  ^  .  f, 

J  sinn  (cu  a)  J 


cu2(cos  cur)  [CsCh2  (cu2a)  ]  dcu 


sinh  (cu  a] 

where  a  =  A2xCn.  We  may  now  integrate  this  integral  with 
respect  to  the  parameter  a. 

J  =  Ida  =  |cu2  cos  (cur)  jcsch2 (cu2a)  dadcu 


‘ 


cos  (cut)  Ctnh  (cu  a)  dcu 


(8.8) 
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We  may  now  use  the  expansion 


Ctnh(aj2a)  =  exp  (A)  +  exp  (-a>2g2 

2  2 
exp  (cd  a)  -  exp  (-(jj  a) 


Thus 


=  (i  +  e~  t0  a)  m  i  +  2  Y]  exp  (-2ka)2u)  (8.9) 

W 


(1  -  e_2ai'a) 


oo 

cos  (cut) 

1  H~ 

2  E 

2 

exp  (-  2ktu  a) 

dm 

J 

GO 

k=l 

_ 

sin  o)T  0 

E 

I  exp 

2 

(-  2km  a)  cos 

asr  da> 

T 

k=l 

J 

s  in  cur  ^ 

oo 

E 

fi(a,7) 

T 

k=l 

(8.10) 


Making  use  of  the  fact  that 


/ 


exp  [-(ax  +  2bx  +  C)]  dx 


1  rr  exp  .  er£  /"Sx  +  b 

2  l  a 


a  *  ^ 

where  the  erf  term  denotes  the  error  function.  We  see  that 


(8.11) 


£2(a>,i)  =  Re  /exp  -  (2ku>  a  -  iair)  dai 


/ 
(i) 


(8.12) 


exp 


8ka 


Re  i  erf 


UrSEaroi  +  — — —  | 

v  2V2ka ^ 


Making  use  of  the  fact  that 

2 


(-l)m  z2m+1 

erf(z)  -  -  /  j  ^ 

\HT  ^  «!  (2m+l ) 


(8.13) 
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and 


^  2ka 


oi  + 


-1 


2yfhka 


=  =  /—  +  2kmi2) 

a  \ 8ka  / 


1/2 


ie 


(8.14) 


where  9  =  tan-  (r/ 2cd)  .  We  obtain  the  results  that 


£1  * 


■5  “  9  o  nt+1/2 

1  ex  _  I__  Vs  (-1)  (Tz/8ka)  +  2km _ (sin  9) 

yj  2ka  8ka  m.'  (2m  +  1) 

(8.15) 


j  sin  cur  _  2 


exp  (-r2/8ko) 

k  =  i  \T^ 


,,  9  9  m+1/2 

(-l)nH‘ir  (t  /8ka)  +  2kato  1 _ sin  9 

m!  (2m  +  1) 


(8.16) 


00  0O 

>-f--2EE 

9a  k-1  m=0 


.m+1 


Cli) - ais_e  exp  (-T28kil> 

m!  (2m  +  1) 


9  9  m+1/2 

f (r  /8ka)  +  2k mi  1 

2ka 
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We  finally  obtain  for  the  correlation  function 
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where  a  =  A  xCn  and  9  =  arc  tan  t/2cu  or  sin  0  =  r/V  4to  +t 

^  I  2  2 

=  t/\  4uj  +t  since  0  is  small.  From  the  definition  of 
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A  =  P/cq2  =  1.2/c0^(  a  =  9.04  (r/f).  From  Eq.  (8.22)  we  obtain 
the  autocorrelation  coefficient  defined  as 


Mil  =  R(r)/R(0) 


In  Ref.  17,  it  was  pointed  out  that  the  total  sound  field  within 
the  Sonic  Fatigue  Facility  when  used  in  the  semi- anechoic  mode 
with  a  reflecting  floor  could  be  predicted  in  terms  of  the  known 
autocorrelation  coefficient  of  the  sound  signal  and  the  radiative 
characteristic  of  the  sound  field.  For  the  case  of  a  simple 
source  radiating  into  infinite  half  space  over  a  hard  reflecting 
surface,  the  expression  giving  the  RMS  pressure  at  any  field 
point  is 


2R(t) 


(8.23) 


In  the  f,ar  field  where  ,  it  is  seen  that  P  depends  upon 

R(t)  and  the  distance  from  the  source. 

For  the  infinitesimal  amplitude  case  this  resulted  in  a 
far  field  with  a  characteristic  interference  pattern.  For  narrow 
bands  of  noise  this  pattern  in  the  direction  of  the  normal  to 
the  floor  varied  from  a  pressure  maximum  at  floor  level  with  suc¬ 
cessive  maximum  and  minimum  pressures  as  height  from  the  floor 
increased.  The  amplitude  of  these  maxima  and  minima  damp  out 
toward  the  ceiling  to  a  level  about  3  db  down  from  that  at  the 
floor.  Discrete  tones  produce  a  similar  pattern  with  the  minima, 
however,  being  very  close  to  zero  pressure  and  the  maxima  6  db 
above  the  level  that  would  exist  at  the  point  for  the  same  source 
radiating  into  free  space.  The  directivity  of  the  source  adds 
sufficient  complexity  to  require  a  computer  program.  If  we  are 
concerned  with  the  near  field,  we  shall  again  require  a  computer 
program  to  use  Eq.  (8.23)  to  predict  the  character  of  the  sound 
field.  It  is  an  assumption  that  Eq.  (8.23)  applies  equally  well 
for  sound  pressure  fields  where  sound  pressure  levels  are  above 
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140  db  re  0.0002  pbar.  If  we  make  this  assumption,  a  computer 
program  could  include  Eq.  (8.23)  and  provide  analytic  predictions 
of  the  sound  field  in  the  Sonic  Fatigue  Facility.  Experimental 
verification  could  be  accomplished  by  sampling  the  sound  field 
in  the  vertical  direction  at  several  locations  in  the  chamber. 
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9. 


EXPERIMENTAL  INVESTIGATION 


The  RTD  siren  operation  will  have  the  capabilities  of 
producing  discrete  frequency  or  narrow  band  noise  signals.  An 
understanding  of  finite  amplitude  propagation  for  both  is,  there¬ 
fore,  desirable.  In  Sections  2-7  we  have  discussed  the  case  of 
discrete  frequency  signals.  In  Section  8,  broad  band  noise  was 
treated  analytically.  In  this  section  experiments  with  band 
limited  noise  signals  are  described. 

9.1  REVIEW  OF  LITERATURE  ON  INTENSE  BROAD  BAND  NOISE 

That  portion  of  the  literature  of  finite  amplitude  sound 
concerned  specifically  with  the  propagation  of  high  intensity 
acoustic  noise  signals  is  quite  limited.  In  our  search  of  the 
literature  we  found  two  papers  by  Ackerman^^  *  .  The  first 

is  an  attempt  at  a  theoretical  treatment  while  the  second  is  a 
limited  attempt  at  the  experimental  verification  of  the  theory. 
The  frequency  band  selected  for  the  experimental  work,  however, 
was  in  the  frequency  range  of  15  kHz  to  50  kHz,  well  above  the 
frequency  range  of  interest  for  the  RTD  Facility.  In  addition, 
overall  sound  pressure  levels  were  in  the  range  of  110  to  135  db, 
while  levels  as  high  as  170  db  are  of  interest  here.  Ackerman's 
theory  predicts  the  effects  of  nonlinearity  of  the  medium  at 
high  sound  pressure  levels  to  be  the  generation  of  sum  and  dif¬ 
ference  frequencies.  Energy  is  thereby  removed  from  the  original 
frequency  band  and  spread  over  larger  parts  of  the  spectrum.  A 
quantitative  description  was  developed  for  the  case  of  1/2  octave 
bands  of  plane  wave  acoustic  noise  three  meters  or  less  from  the 
generator.  At  15,000  Hz  this  distance,  three  meters,  corresponds 
to  130  wavelengths. 

The  main  siren  bank  at  the  RTD  Facility  produces  acous¬ 
tic  energy  of  significant  amount  in  the  frequency  range  of 
approximately  150  to  2000  Hz,  In  this  frequency  range  sound 
pressure  levels  as  high  as  165  db  re  0.0002  iibar  for  a  bandwidth 
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of  50  Hz  have  been  measured  in  the  near  field  of  the  sirens. 

For  frequencies  in  the  lower  part  of  this  range,  150  to  500  Hz 
and  at  levels  below  170  db,  essentially  two  regimes  are  possible, 
according  to  the  theory.  At  sound  pressure  levels  of  up  to 
130  db,  infinitesimal  solutions  are  valid,  that  is,  no  signi¬ 
ficant  finite  amplitude  effects  are  expected  to  occur.  At 
sound  pressure  levels  above  130  db  "sum  frequencies"  build  up 
but  a  shock  wave  does  not  form  in  three  meters  even  at  levels 
of  up  to  170  db.  In  the  second  frequency  range,  500  Hz  to 
2000  Hz,  the  upper  limit  at  which  infinitesimal  theory  is  still 
valid  is  lowered  to  about  110  to  120  db  SPL.  The  lower  level 
occurs  at  the  higher  frequency.  At  sound  pressure  levels  be¬ 
tween  these  levels  and  about  155  db  "sum  frequencies"  build  up 
although  no  shock  wave  is  possible.  Above  this  level  the  theory 
fails  to  offer  a  prediction. 

The  high  frequency  sirens  can  produce  significant  acous¬ 
tic  energy  in  the  frequency  range  between  500  and  9600  Hz.  It 
is,  therefore,  of  interest  to  consider  the  nonlinear  propaga¬ 
tion  of  acoustic  signals  up  to  9600  Hz  as  well.  For  frequen¬ 
cies  above  2000  Hz,  the  same  regimes  as  those  from  500  Hz  to 
2000  Hz  also  occur.  However,  the  sound  pressure  level  regime 
for  each  successive  frequency  range  is  correspondingly  lower. 
Thus,  one  is  to  expect  the  growth  of  "sum  frequencies"  at  sound 
pressure  levels  as  low'as  105  db.  The  theory  again  fails  at 
levels  above  155  db  at  these  frequencies.  The  theory  adds  one 
further  prediction  for  noise  bands  at  these  higher  frequencies. 

At  sound  pressure  levels  of  150  db  shock  formation  cannot  occur 
because  of  excessive  attenuation  of  the  "sum  frequencies"  gen¬ 
erated  by  nonlinearities  in  the  medium.  It  is  to  be  noted  that 
the  theory  offers  no  prediction  at  levels  above  150  db  in  the 
upper  end  of  this  frequency  range,  about  500  Hz  to  10,000  Hz. 

The  experimental  results  have  been  confined  to  the  high 
sonic  and  low  ultrasonic  frequency  range.  In  this  range,  the 
excessive  absorption  of  the  medium  for  the  fundamental  and 
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whatever  "sum  frequencies"  are  generated  is  the  major  factor. 

As  a  consequence,  shock  formation  is  precluded  in  the  theory 
even  at  very  high  sound  pressure  levels  up  to  170  db.  At  the 
lower  sonic  frequencies  of  interest  here,  significant  energy 
can  be  anticipated  in  the  "sum  frequencies"  and  the  propagation 
of  bands  of  noise  at  the  higher  sound  pressure  levels  can  be 
expected  to  be  significantly  different. 

9.2  EXPERIMENTAL  PROCEDURES 

The  nonlinearity  of  the  basic  equations  applicable  to 
sound  propagation  at  these  sound  pressure  levels  implies  that 
the  simple  summation  procedures  of  the  linear  theory  are  not 
per  se  applicable.  Thus,  one  cannot  say  a  priori,  that  integra¬ 
tion  of  pure  tones  over  the  band  of  interest  constitutes  a  solu¬ 
tion  to  the  problem  in  the  nonlinear  case.  It  is  for  these 
reasons  that  we  have  undertaken  an  experimental  investigation 
of  the  propagation  of  plane  waves  of  acoustic  noise. 

The  major  experimental  limitation  was  the  generation  of 
noise  signals  at  the  sound  pressure  levels  of  interest  (up  to 
170  db)  without  initially  introducing  'a  significant  level  of  har¬ 
monic  distortion.  This  limitation  required  that  the  sound  be 
generated  with  multiple  transducers  each  operating  at  a  low  power 
level.  The  low  power  level  per  transducer  required  a  large  num¬ 
ber  of  transducers  whose  sound  energy  output  is  converged  to  a 
narrow  channel  to  raise  the  intensity  to  the  degree  required. 

For  this  purpose  use  was  made  of  the  high  intensity  pro¬ 
gressive  sound  tube  facilities  at  the  Riverbank  Acoustical 
Laboratory  in  Geneva,  Illinois.  These  facilities  are  described 
in  Ref.  20  and  briefly  summarized  in  Appendix  B.  For  this  pro¬ 
gram  the  low  frequency  and  medium  frequency  impedance  tubes  were 
used . 

The  low  frequency  tube  used  has  a  usable  length  of  20  ft 
(6.10  meters)  and  a  diameter  of  3  in.  (7.5  cm).  This  was  terminated 


55 


by  a  non-reflecting  termination  consisting  of  a  wedge  of  acous¬ 
tically  absorbing  material.  The  cut  off  frequency,  f c ,  the 
frequency  of  the  first  cross  mode,  for  a  cylindrical  tube  of 
radius,  a,  is  given  by  Morse as: 

f  =  0.586  c/2a 

For  a  3  in.  tube  the  cut  off  frequency  is  2650  Hz.  No  difficulty 
was  experienced  in  propagating  bands  of  noise  down  the  tube  at 
levels  up  to  150  db  SPL  at  frequencies  as  high  as  1000  Hz. 

The  sound  source  for  this  tube  is  a  bank  of  49  speakers 
arranged  in  a  square  array  of  7  by  7,  The  speakers  were  8  in. 
cone  loudspeakers,  Jensen  P8-PC.  The  arrangement  of  the  speakers, 
their  power  handling  capabilities,  and  efficiency  are  discussed 
in  Ref.  20  together  with  construction  and  acoustical  performance 
of  the  exponential  horn  coupling  the  speakers  to  the  tube. 

In  order  to  attain  higher  levels  at  frequencies  from 
100  Hz  to  2000  Hz  we  used  the  medium  frequency  tube.  The  de¬ 
tails  of  this  facility  are  also  discussed  in  Ref.  20.  The  modi¬ 
fication  to  this  facility  made  for  this  experiment  included: 

(a)  adding  a  pc  termination,  (b)  adding  a  3  ft  section  to  the 
working  length  of  the  tube,  and  (c)  replacing  the  power  ampli¬ 
fier  with  the  Savage  1  kilowatt  audio  amplifier.  With  10  per¬ 
cent  band  limited  noise  having  center  frequencies  1500  Hz  and 
2000  Hz  could  be  produced  at  sound  pressure  levels  of  130  db 
re  0.0002  pbar. 

In  each  tube  the  microphone  can  be  moved  down  the  axis 
of  the  tube  with  its  diaphragm  normal  to  this  axis.  Our  pro¬ 
cedure  was  to  excite  the  drivers  with  a  10  percent  band  limited 
noise,  position  the  microphone  at  a  known  distance  from  the 
cross  section  where  the  pc  ending  began,  and  to  record  the 
microphone  output  on  a  Magnecorder,  Model  PT6BA  magnetic  tape 
recorder.  A  calibration  signal  of  a  pure  sinusoid  at  a  known 
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SPL  and  at  the  center  frequency  of  the  noise  band  was  also 
recorded  on  the  tape . 

Simultaneously  with  the  recording  of  the  noise  signal 
at  a  given  station  this  station  was  identified  on  a  second 
channel.  These  tapes  were  made  into  tape  loops,  one  for  each 
station  at  which  analysis  was  to  be  performed.  The  loops  were 
played  through  the  Ampex  recorder  and  the  sound  pressure  level 
in  a  10  Hz  band  relative  to  that  of  the  calibrate  signal  was 
measured  using  the  Quan-Tech  Laboratories  Model  303  wave  analyzer 
This  was  done  in  two  ways.  In  the  first,  each  harmonic  of  the 
fundamental  on  a  given  tape  loop  was  measured  using  the  meter 
reading  of  the  wave  analyzer  directly.  The  peak  level  in  the 
frequency  region  around  the  harmonic  was  sought  and  recorded. 

In  a  second  approach  the  tape  recorder  signal  was  pre¬ 
sented  to  the  Quan-Tech  wave  analyzer  whose  output  was  presented 
to  a  Bruel  and  Kjaer  Type  2305  level  recorder.  The  level  re¬ 
corder  was  driven  synchronously  by  the  Quan-Tech  wave  analyzer. 
Writing  speed  and  paper  speed  on  the  level  recorder  were  adjusted 
to  provide  maximum  detail  compatible  with  stability  and  repro¬ 
ducibility.  A  five  second  tape  loop  was  completely  analyzed 
over  a  frequency  range  of  from  200  and  10,000  Hz  in  approximately 
one  hour.  When  empty  reels  of  tape  were  played  through  this 
system,  the  noise  level  in  a  10  Hz  band  was  at  least  70  db  down 

from  the  peak  level  of  the  fundamental  at  all  frequencies  above 
250  Hz. 

It  is  desirable  to  have  a  signal  output  from  the  drivers 
and  the  exponential  horn  as  free  as  possible  from  harmonic  dis¬ 
tortion.  To  establish  how  well  this  requirement  had  been  achieved 
we  probed  the  sound  field  immediately  in  front  of  the  speaker 
bank  and  for  a  distance  of  42  in.  into  the  horn  from  the  plane 
of  the  speakers.  A  10  percent  bandwidth  of  noise  with  a  center 
frequency  of  250  Hz  at  a  level  of  157  db  was  generated. 
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Sound  pressures  were  recorded  as  a  function  of  distance 
from  the  plane  of  the  speakers  along  the  axis  of  the  horn. 
For  each  of  eight  microphone  positions,  the  recordings 
were  made  into  tape  loops  and  analyzed.  The  second  and 
third  harmonic  distortion,  that  is,  the  difference  in  level 
between  the  first  and  second  and  between  first  and  third 
harmonic  components  at  each  point  expressed  in  db  is 
plotted  as  a  function  of  distance  in  Fig.  9.1.  Distortion 
products  follow  a  standing  wave  pattern  similar  to  the 
standing  wave  of  the  fundamental  (not  shown). 

We  assume  that  strong  reflections  occur  at  the 
transition  from  the  horn  throat  to  the  constant  diameter 
experimental  section.  The  cross-section  is  continuous  in 
this  transition  region,  but  the  rate  of  change  of  section 
is  discontinuous.  Consequently  there  is  an  impedance  dis¬ 
continuity  giving  rise  to  reflections  and  thereby  standing 
waves.  The  variation  in  amplitude  of  the  fundamental,  as 
well  as  that  of  the  second  and  third  harmonic  follow  the 
typical  standing  wave  pattern  usually  observed.  The  wave¬ 
length  corresponds  to  the  fundamental  frequency  with  only 
second  order  perturbations  at  the  harmonic  frequency.  More¬ 
over  the  harmonic  distortion  appears  to  predominate  at  the 
pressure  maxima.  This  leads  to  the  supposition  that  the 
harmonics  are  mainly  generated  in  the  concentrator  horn 
rather  than  in  the  transducers  or  amplifiers.  A  reduction 
of  standing  wave  patterns  in  the  coupling  horn  might  be 
achieved  by  changes  in  the  horn  design  and  would  lead  to 
some  slight  reduction  in  the  initial  distortion  in  the 
measuring  section.  Such  modifications  were  outside  the 
scope  of  the  present  project  and  some  initial  distortion 
in  the  test  section  had  to  be  accepted. 
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Fig.  9.1.  Level  of  Second  and  Third  Harmonic  Components 
below  the  Level  of  the  Fundamental  as  a 
Function  of  Distance  from  the  Speakers  in  the 
Horn  Coupling  Section.  (Band  limited  noise 
bandwidth  10  percent  mean  frequency  250  Hz 
at  157  db  re  0.0002  pbar  near  speakers.) 
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The  measurements  made  to  establish  the  extent  of 
harmonic  distortion  introduced  by  the  speakers  and  by  non¬ 
linear  propagation  effects  in  the  horn  were  not  conclusive. 
However,  harmonic  distortion  in  the  signal  at  the  entrance 
to  the  acoustic  tube  was  sufficiently  low  so  that  the  growth 
of  harmonics  could  be  observed  and  measurements  could  be 
made. 

Measurements  of  harmonic  content  of  the  signal  in 
the  tube  entrance  at  the  horn  transition  for  a  band  of  noise 
with  central  frequency  at  800  Hz  and  overall  SPL  of  145  db 
showed  the  second  harmonic  down  30  db  from  the  fundamental 
and  the  third  harmonic  down  50  db.  Levels  of  fourth  and 
higher  harmonics  were  below  the  noise  threshold  of  available 
ins  trumenta t ion . 

The  lowest  central  frequency  used  for  the  narrow 
band  (,10  percent  bandwidth)  noise  signals  in  the  study  was 
300  Hz.  Figures  9.2  and  9.3  relate  to  the  measurements 
made  at  this  frequency.  Figure  9.2  presents  spectral  analy¬ 
sis  of  the  recordings  of  the  noise  signal  at  a  distance  of 
115  in.  from  the  entrance  to  the  tube.  The  overall  noise 
level  for  Fig.  9.2  was  130  db.  For  a  level  of  110  db  no 
distortion  was  measurable.  In  Fig.  9.3  we  have  presented 
the  analysis  at  the  entrance  to  the  tube,  at  the  distance 
of  115  in.  and  at  a  distance  of  230  in.  The  overall  sound 
pressure  level  was  160  db. 
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9.2.  Spectrum  of  Sound  of  a  10  percent  Band  of 
Noise  Analyzed  with  a  10  Hz  Bandwidth 
Analyzer.  Microphone  115  in.  from  start 
of  straight  section  of  acoustic  tube. 


61 


62 


At  the  overall  level  of  130  db,  Fig.  9.2  ,  energy  is 
present  in  the  second  and  third  harmonic  but  their  levels  are 
35  db  and  42  db  lower  than  the  fundamental.  At  an  overall 
level  of  160  db  (Fig.  9.3),  however,  there  is  measurable  energy 
at  even  the  seventh  harmonic,  while  the  lower  harmonics,  second 
and  third,  have  levels  within  8  db  and  19  db  respectively  of 
the  fundamental  at  a  distance  of  115  in.  from  the  speakers. 

The  high  level  of  energy  in  the  second  and  third  harmonics  at 
the  position  nearest  the  opening  to  the  tube  from  the  horn  to¬ 
gether  with  significant  energy  in  the  fourth  and  fifth  harmonics 
is  not  understood.  Such  results  were  not  observed  at  higher 
frequencies  at  the  very  highest  overall  levels.  The  investiga¬ 
tion  of  the  sound  field  in  the  exponential  horn  at  comparable 
overall  levels,  157  db,  and  at  a  frequency  of  250  Hz  did  not 
suggest  this  result.  Even  this  degree  of  harmonic  content,  how¬ 
ever,  represents  only  a  small  distortion  of  the  noise  signal. 

It  cannot  be  said,  for  example,  that  a  shock  wave  is  present. 

In  Fig.  9.4  we  have  plotted  the  variation  of  sound  pres¬ 
sure  level  in  all  harmonics  present  versus  distance  down  the 
tube  from  the  entrance.  The  central  frequency  is  300  Hz  and 
the  overall  level  is  160  db.  At  this  frequency  the  working 
length  of  the  tube  was  about  1.6  wavelengths.  There  appears  to 
be  no  significant  change  in  level  with  distance  for  this  dis¬ 
tance  in  the  first  four  harmonics.  The  fifth  and  sixth  harmonics 
have  increased  in  level  about  2  to  3  db  over  1.6  wavelengths. 

The  appearance  and  growth  to  a  maximum  level  of  harmonics 
in  a  noise  signal  is  well  illustrated  in  Fig.  9.5.  Here  the 
central  frequency  was  800  Hz  and  the  overall  level  for  the  sig¬ 
nal  was  145  db.  This  was  the  maximum  level  attainable  in  the 
tube  at  this  frequency.  The  wavelength  of  sound  at  800  Hz  is 
0.43  meters  so  that  the  working  lengths  of  the  tube  correspond 
to  13.7  wavelengths.  The  appearance  of  successively  higher 
harmonics  and  their  rapid  growth  to  peak  levels  as  the  wave 
progresses  down  the  tube  is  a  prominent  feature.  The  levels  of 
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Growth  of  Harmonics  of  a  Plane  Wave  of  Narrow  Band  Noise  in  an 
Impedance  Tube  (Overall  sound  pressure  level  in  10  percent 
bandwidth  160  db,  central  frequency  300  Hz). 
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Growth  of  Harmonics  of  a  Plane  Wave  of  Narrow  Band  Noise  in  an 
Impedance  Tube  (Overall  sound  pressure  level  in  10  percent 
bandwidth  145  db,  central  frequency  800  Hz). 


the  second  and  third  harmonics  at  the  entrance  to  the  tube  from 
the  exponential  horn  are  30  db  and  50  db  down  from  the  funda¬ 
mental,  Therefore,  the  possible  distortion  produced  in  the  sig¬ 
nal  by  the  speakers  and  the  horn  do  not  obscure  the  effects  of 
nonlinear  propagation  in  the  tube.  At  overall  levels  of  this 
order  the  levels  attained  by  harmonics  higher  than  the  third 
will  not  be  significant  for  operation  of  the  Sonic  Fatigue 
Facility.  Within  the  tube  these  higher  harmonics  rise  rapidly 
initially  but  over  the  distance  of  four  to  five  wavelengths  have 
already  passed  their  peak.  The  maximum  level  attained  for  the 
second  harmonic  was  16  db  below  the  fundamental  and  for  the 
third  harmonic  was  30  db  below  the  fundamental. 

The  uniformity  of  the  sound  pressure  level  of  the  funda¬ 
mental  along  the  length  of  the  tube  attests  to  the  quality  of 
the  termination.  The  standing  wave  ratio  seems  to  have  been 
below  3  db.  The  quarter  wavelength  peak  to  valley  distance  is 
approximately  0.1m  for  800  Hz,  hence  the  apparent  irregularity 
of  the  point  shown  in  Fig.  9.5,  some  of  which  come  close  to 
maxima  and  some  are  near  minima  of  the  small  residual  standing 
wave  pattern.  This  residual  standing  wave  pattern  tends  to 
obscure  the  growth  of  harmonics  near  the  far  end  of  the  tube. 
Indications  are  strong  that  the  rate  of  growth  of  the  harmonics 
decreases  in  the  manner  indicated  for  pure  tones  in  Fig.  2.3. 
Comparison  with  pure  tone  assumptions  is  shown  dotted  in  Fig. 

9,5  and  discussed  together  with  a  similar  curve  at  1000  Hz  in 
Section  9.3. 

In  plane  wave  propagation  the  power  of  harmonics  is 
less  than  2  percent  of  the  total  after  travelling  a  distance  of 
fourteen  wavelengths  at  an  initial  fundamental  sound  pressure 
level  of  145  db  at  800  Hz.  This  two  percent  reduction  in  funda¬ 
mental  level  is  not  discernible  among  the  residual  standing 
wave  fluctuations . 

These  results  pertain  to  a  signal  confined  to  the 
interior  of  a  tube.  They,  therefore,  represent  upper  bounds  for 
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the  case  of  the  sirens  of  the  Sonic  Fatigue  Facility  where  spheri 
cal  spreading  will  decrease  intensity  within  the  wave  as  it  props 
gates  out  into  the  room.  With  spherical  spreading  propagations 
of  the  wave  for  a  distance  of  6  meters  would  decrease  the  sound 
pressure  level  in  the  signal  14  db  below  its  level  at  1  meter 
from  the  source.  This  applies  to  all  harmonics  as  well.  It  is 
directional ,  nevertheless,  spreading  of  the  wave  energy  must 
occur  and  the  transfer  of  energy  into  the  higher  harmonics  must 
be  less  efficient  than  is  true  for  the  acoustic  tube. 

Figure  9.6  presents  comparable  results  for  a  band  of 
noise  10  percent  bandwidth  and  centered  at  1000  Hz.  Again  the 
second  harmonic  has  reached  a  peak  at  a  level  16  db  below  the 

fundamental  and  the  third  harmonic  has  reached  a  peak  at  a  level 

26  db  below  the  fundamental. 

At  sound  pressure  levels  of  130  db  re  0.0002  p.bar  or 
lower  for  frequencies  below  500  Hz  and  at  110  db  or  lower  for 

the  frequencies  of  800  Hz  and  1000  Hz,  there  was  no  evidence 

of  a  growth  of  harmonic  content  as  the  noise  signal  propagated 
down  the  tube.  For  these  frequencies  and  at  these  levels  the 
infinitesimal  theory  clearly  holds.  Although  a  small  harmonic 
content  is  present  in  the  s ignal  at  the  entrance  to  the  tube , 
it  does  not  increase  at  any  harmonic.  These  results  are  in 
agreement  for  these  frequencies  with  Ackerman's  theory.  We  have 
presented  a  resume  of  the  data  for  these  cases  in  Table  9.1, 

Data  for  three  stations  are  included  (a)  near  the  termination, 

(b)  near  the  half  way  point  in  the  tube,  and  (c)  at  the  speaker 
end  of  the  tube. 

At  the  maximum  levels  achieved  in  the  tube  and  the  path 
length  of  travel  for  the  signal,  no  shock  wave  could  be  said 
to  have  formed.  This  is  true  even  for  the  case  of  the  signal 
with  central  frequency  of  300  Hz  and  overall  sound  pressure 
level  of  160  db.  There  was  a  significant  transfer  of  energy 
into  the  second  and  third  harmonics  at  the  higher  overall  sound 
pressure  levels  above  145  db  which  could  be  attributed  to 
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Fig.  9.6.  Growth  of  Harmonics  of  a  Plane  Wave  of  Narrow  Band  Noise  in  an 
Impedance  Tube  (Overall  sound  pressure  level  in  10  percent 
bandwidth  145  db,  central  frequency  1,000  Hz). 


nonlinear  effects  in  propagation  of  the  signal  through  the  tube. 
The  maximum  transfer  of  energy  into  these  harmonics,  however, 
brought  them  to  levels  no  higher  than  2  percent  for  the  second 
harmonic  and  0.1  percent  for  the  third  harmonic  of  the  energy 
in  the  fundamental.  For  the  higher  harmonics  the  energy  con¬ 
tent  at  the  most  is  negligible.  The  harmonic  distortion  in 
the  sound  output  of  the  sirens  will  be  expected  to  be  much  more 
significant  than  the  additional  transfer  of  energy  into  the 
higher  harmonics  attributable  to  nonlinear  propagation  in  the 
volume  of  the  room. 


TABLE  9.1 

LIMITED  GROWTH  OF  HARMONICS  IN  10  PERCENT  BAND  LIMITED  NOISE 
FOR  RELATIVELY  LOW  SOUND  PRESSURE  LEVELS 


Distance  from 
Start  of 


db  Below  Fundamental 


’requency 

(Hz) 

SPL  db 

re  0.0002pbar 

Straight  Section 
(inch) 

f-h2 

f-h3 

F-H4  f-h5 

f-h6 

300 

110 

0 

32 

37 

39  41 

43 

120 

30 

34 

36  38 

39 

236 

28 

31 

35  37 

38 

300 

130 

0 

37 

44 

54  58 

65 

120 

36 

43 

54  58 

63 

236 

37 

45 

54  58 

61 

500 

110 

0 

34 

45 

49 

120 

34 

43 

50 

236 

36 

44 

50 

500 

130 

0 

32 

47 

60 

120 

32.5 

47.5 

58.5 

236 

29 

47 

58 

800 

110 

0 

35.5 

46.5 

120 

37 

49 

236 

37 

51 

1000 

100 

0 

43 

46 

52 

120 

44 

48 

54 

236 

44 

48 

53 
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9.3  COMPARISON  OF  EXPERIMENTAL  RESULTS  WITH  COMPUTATIONS 


It  was  pointed  out  in  Section  2,6  that  distortion  effects 
may  be  expected  in  band  limited  noise  in  the  same  way  as  in  pure 
tones,  but  the  distortion  may  bear  a  different  relationship  to 
the  root  mean  square  amplitude  because  of  differences  in  the  dis¬ 
tribution  of  peak  amplitudes.  In  Fig.  9.5  the  growth  of  the 
second  and  third  harmonic  of  a  pure  tone  of  the  same  RMS  ampli¬ 
tude  is  indicated  and  seen  to  be  very  close  to  the  amplitude  of 
the  distortion  products  measured  for  the  noise.  Such  agreement 
might  be  fortuitous.  In  Fig.  9.6  the  growth  of  harmonics  of  a 
pure  tone  7  db  higher  in  intensity  than  the  noise  is  shown.  The 
measured  rate  of  growth  is  seen  to  be  less  than  that  computed 
for  the  tone.  If  a  lower  level  tone  had  been  chosen  for  com¬ 
parison  a  lower  growth  rate,  but  also  a  lower  level  of  distor¬ 
tion  products  would  have  been  obtained.  It  is  possible  that  the 
discrepancy  may  be  accounted  for  in  part  by  the  finite  initial 
distortion. 

Further  theoretical  and  experimental  investigations 
along  this  line  would  seem  appropriate. 
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10.  ROOM  TREATMENT 


It  was  concluded  in  Sections  2-9  that  effects  in  the 
room  at  points  remote  from  the  siren  are  relatively  small  com¬ 
pared  with  nonlinear  effects  in  the  siren  horns. 

The  interactions  due  to  reflections  at  the  floor  for 
the  semi-anechoic  mode  of  operation  and  due  to  the  reflections 
at  walls  and  floor  in  the  semi- reverberant  mode  may  be  computed 
from  the  treatment  in  Section  7.  For  normal  incidence  where 
long  wave  trains  are  superimposed  over  considerable  distances 
slight  additional  distortion  and  intermodulation  is  predicted. 
However,  the  shape  of  the  room  and  the  absence  of  a  wall  paral¬ 
lel  to  the  plane  of  the  sirens  limits  normal  incidence  to  an 
insignificant  portion  of  the  total  sound  energy.  For  oblique 
incidence  the  interaction  will  be  only  over  a  limited  region, 
where  incident  and  reflected  sound  waves  overlap.  Moreover, 
positive  and  negative  pressure  regions  overlap  alternately  so 
that  the  overall  effect  tends  to  cancel  out.  This  means  that 
the  nonlinear  effects  are  simply  those  of  direct  propagation 
from  an  effective  image  source.  The  contribution  of  reflected 
waves  to  nonlinear  distortion  and  intermodulation  is  negligible 
for  practical  purposes;  because,  as  previously  stated,  nonlinear 
propagation  effects  in  the  spherical  propagation  region  are 
small  in  the  frequency  and  intensity  regions  under  considera¬ 
tion,  and  also  because  the  first  reflections  are  all  near  graz¬ 
ing  incidence. 

Thus,  room  conditions  whether  reverberant  or  anechoic  do 
not  change  any  of  the  basic  contributions  of  nonlinear  effects. 
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11.  RECOMMENDATIONS 


11.1  FACILITY  USE 

The  harmonic  content  of  the  sound  should  be  taken  into 
account  when  setting  up  specific  sound  spectra  in  the  room. 
Direct  generation  of  high  frequency  components  may  be  reduced 
due  to  the  generation  of  these  components  by  distortion  in  the 
lower  frequency  component  generators . 

11.2  MEASUREMENTS  AND  TESTS 

The  mathematical  treatment  predicts  the  level  of  the 
harmonic  components  for  finite  amplitude  generators.  However, 
at  the  highest  amplitudes  the  extrapolation  is  somewhat  uncer¬ 
tain  and  direct  measurements  are  recommended.  Some  of  these 
measurements  have  already  been  undertaken.  It  is  recommended 
that  these  measurements  be  extended  to  a  more  detailed  mapping 
of  the  sound  pattern  of  a  single  siren.  The  amplitude  of  both 
fundamental  and  of  each  harmonic  should  be  measured  at  a  con¬ 
venient  distance  from  the  source  as  a  function  of  the  angle  be¬ 
tween  the  microphone  direction  and  the  horn  axis.  Two  planes 
should  be  chosen,  preferably  the  horizontal  plane  and  the  plane 
containing  the  diagonal  of  the  square  horn  mouth.  If  the  direc¬ 
tional  characteristics  of  the  sirens  is  found  to  be  very  differ¬ 
ent  in  these  two  planes  so  that  interpolation  is  uncertain, 
intermediate  planes  inclined  at  22.5°,  11.25°,  and  37.75°  should 
be  chosen  to  augment  the  data. 

These  data  will  be  needed  if  precise  fields  are  to  be 
computed  in  the  facility.  The  directional  characteristics  of 
the  sirens  may  be  amenable  to  computation.  In  that  case  a  more 
limited  experimental  program  can  be  set  up  in  which  these 
characteristics  are  measured  at  a  very  limited  number  of  fre¬ 
quencies  to  check  the  computations. 
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11.3 


FACILITY  MODIFICATIONS 


11.3.1  Low  Frequency  Horns 

Facility  modifications  have  been  contemplated  to  extend 
the  capability  to  lower  frequencies  by  modifying  some  of  the 
horns  to  have  a  lower  cut  off  frequency.  This  modification 
should  be  limited  to  the  minimum  number  of  horns  to  be  used  at 
low  frequency.  At  small  amplitudes  the  cut  off  frequency  has 
negligible  effect  on  high  frequency  performance.  At  large  am¬ 
plitudes  the  distortion  is  greater  if  the  sound  is  maintained 
at  the  high  intensity  for  a  long  time  in  a  longer  horn. 

11.3.2  High  Frequency  Horns 

Some  consideration  may  also  be  given  to  the  relative 
merits  of  reducing  the  high  frequency  distortion  by  Increasing 
the  flare  in  some  of  the  horns.  The  cut  off  frequency  will 
thereby  be  increased,  thus  adversely  affecting  the  flexibility 
of  the  facility.  A  compromise  must  be  established  taking  into 
account  the  expected  use  of  the  facility.  In  other  words,  if 
the  distortion  is  considered  unimportant,  and  a  need  for  high 
intensity  at  low  frequency  is  foreseen,  a  larger  number  of  low 
frequency  units  would  be  .desirable ,  than  if  maximum  output  with 
low  distortion  is  required  at  high  frequencies. 

11.4  FURTHER  THEORETICAL  ANALYSIS 

The  importance  of  the  distortion  produced  in  the  horn 
of  the  siren  (Section  3)  was  not  realized  sufficiently  early 
in  the  program  to  devote  the  time  necessary  to  extend  existing 
theories  to  the  very  large  amplitudes  involved.  Existing 
theories  were  cautious ly  extrapolated ,  but  the  precise  computa¬ 
tion  of  the  value  of  the  distortion  products  requires  a  more 
sophisticated  theory.  The  theory  used  is  essentially  a  second 
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order  approximation  and  is  valid  for  moderate  distortion  levels. 
The  theory  needs  to  be  taken  to  a  higher  order  of  approximation 
to  establish  the  precise  level  of  distortion  products.  This 
knowledge  would  permit  a  more  exact  evaluation  of  the  relative 
merits  of  using  horns  specifically  designed  for  narrow  frequency 
regions  as  discussed  under  11.3. 

11.5  EXPERIMENTAL  STUDIES 

11.5.1  Distortion  in  Siren  Horns 

The  portions  of  the  curves  shown  dotted  in  Fig.  3.1  may 
also  be  verified  or  modified  by  an  experimental  study.  A  simple 
series  of  tests  would  consist  of  placing  a  suitable  transducer 
close  to  the  mouth  of  a  siren  and  measuring  the  amplitude  of 
each  of  the  harmonics  generated  (including  the  fundamental)  as 
a  function  of  power  input  to  the  siren.  The  experiment  would 
be  repeated  for  different  frequencies  over  the  full  range  of  the 
siren  capability.  A  study  of  this  kind  is  relatively  easy  to 
undertake  and  is  strongly  recommended  to  augment  or  replace  the 
theoretical  study  of  Section  11.2. 

In  addition  it  would  be  useful  if  the  directional  pat¬ 
terns  of  the  radiation  of  the  sound  from  the  sirens  be  deter¬ 
mined  experimentally.  This  will  be  needed  if  the  sound  field 
in  the  facility  is  to  be  computed.  The  question  of  the  radiation 
of  a  sound  of  a  certain  frequency  may  be  affected  by  the  pres¬ 
ence  of  other  sounds.  Thus  the  radiation  of  a  3  kHz  component 
when  generated  as  a  fundamental  may  be  different  from  the  radia¬ 
tion  pattern  of  the  sound  of  this  frequency  when  it  is  the  second 
harmonic  of  1.5  kHz  or  the  third  harmonic  of  1  kHz. 

11.5.2  Other  Nonlinear  Effects 

An  experimental  investigation  of  other  effects  studied 
in  this  report  may  also  be  of  interest.  This  includes  the 
reflection  of  simple  large  amplitude  sounds  and  the  interaction 
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of  sounds  of  different  frequency  generating  intermodulation 
products  at  the  sum  and  difference  frequencies.  The  study  of 
the  precise  level  of  distortion  in  horns  (Section  11.5.1)  should 
be  given  high  priority. 

11.5.3  Noise 

The  magnitude  of  the  nonlinear  effects  on  initially 
simple  harmonic  pressure  fluctuations  is  not  readily  translated 
to  random  signals;  Sections  2.5,  8,  and  9.3  indicate  some  of 
the  difficulties  involved.  It  is  suggested  that  an  experimental 
program  be  set  up  in  the  high  intensity  facility  in  which  ampli¬ 
tudes  of  the  fundamental  and  the  first  few  harmonics  be  measured 
at  various  distances  on  the  axis  of  a  siren  at  different  levels 
of  pure  tones  and  narrow  band  noise.  Simultaneous  theoretical 
development  as  in  Section  8  would  also  be  appropriate. 
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12.  CONCLUSIONS 


The  nonlinearities  introduced  by  the  medium  for  finite 
amplitude  sound  waves  introduce  harmonics  into  the  sound  waves 
generated  by  the  sirens.  This  means  that  in  order  to  generate 
a  predetermined  spectrum  in  the  facility  allowance  must  be 
made  for  this  effect.  The  effect  can  be  reduced  by  redesigning 
some  of  the  horns  to  have  higher  cut  off  frequencies.  In  any 
case  if  lower  cut  off  frequency  horns  are  to  be  used  on  some 
of  the  sirens  in  order  to  generate  lower  frequency  sounds,  the 
changeover  should  be  limited  to  the  number  needed  for  the  low 
frequency  generation.  The  use  of  a  horn  with  a  low  cut  off 
frequency  for  high  frequency  generation  will  increase  the  dis¬ 
tortion  produced. 

Because  of  the  large  unavoidable  distortion  in  the  siren 
horns,  the  much  smaller  additional  effects  in  the  room  may  be 
neglected. 
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APPENDIX  A 


MEASUREMENTS  OF  SOUND  FIELDS  IN  THE  AFFDL  FACILITY 

Sound  pressure  levels  were  measured  in  the  Facility  with 
siren  #2  only,  siren  #4  only,  and  both  siren  #2  and  siren  #4 
operating  together.  Waveforms  were  observed  for  operating  fre- 
quenc ies  of  300  Hz  and  1000  Hz .  At  300  Hz  distortion  was  re¬ 
ported  to  be  low  and  detailed  measurements  were  not  further  re¬ 
ported. 

For  a  frequency  of  1000  Hz  waveforms  of  the  sound  were 
recorded  immediately  in  front  of  siren  #2  and  at  several  loca¬ 
tions  in  the  room.  The  locations  were  chosen  in  a  single  plane 
normal  to  the  plane  of  the  sirens  midway  between  #2  and  #4  at 
distances  of  6  ft.  2  in.,  12  ft.  8.5  in.,  and  15  ft.  4  in.  and 
heights  of  3  ft,  11  in.,  10  ft.,  12  ft,,  14  ft.,  and  16  ft.  above 
the  floor. 

Spectrum  analysis  yielded  a  rather  complex  picture  of 
the  distortion  at  various  points.  This  is  attributed  to  the 
fact  that  the  radiation  pattern  of  the  sirens  is  rather  complex 
and  the  angle  of  the  siren  axis  and  the  vector  from  the  mouth 
of  the  siren  to  the  measuring  location  changes  with  the  measur¬ 
ing  location,  so  that  peaks  or  valleys  of  the  radiation  pattern 
may  be  close  to  different  measuring  points  for  the  different 
harmonics . 

The  sound  pressure  level  in  front  of  siren  #2  was  168  db 
and  the  2nd,  3rd,  and  4th  harmonic  were  10  db,  14  db,  and  15  db 
respectively  below  the  fundamental. 

The  largest  distortion  in  the  room  was  measured  at  the 
12  ft.  height  above  the  floor  at  the  12  ft.  8.5  in.  distance 
where  the  fundamental  was  150  db  and  the  2nd,  3rd,  and  4th  har¬ 
monic  were  only  5  db,  10  db,  and  11  db  respectively  below  the 
fundamental;  at  most  other  locations  the  distortion  was  no  greater 
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than  immediately  in  front  of  the  horn.  It  seems  reasonable  to 
explain  this  behavior  by  assuming  that  the  harmonics  showed 
side  lobes  in  their  radiation  pattern  in  this  direction. 

The  recorded  waveforms  are  shown  in  Fig.  A-L,  Very  dis¬ 
tinct  triangular  shock  waves  are  readily  discerned.  Some  of  the 
records  have  been  inverted  because  no  attention  had  been  given 
to  the  phase  of  the  microphone  connection  and  the  presentation 
of  the  shock  with  uniform  polarity  helps  distinguish  this  fea¬ 
ture  in  the  record.  The  columns  of  oscillograms  show  simultaneous 
records  of  six  microphones  at  a  distance  from  the  siren  bank 
shown  at  the  top  of  the  column.  The  height  of  the  microphone 
above  the  floor  is  shown  in  the  left  margin.  The  fourth  oscillo¬ 
gram  is  the  record  of  a  monitor  microphone  located  in  the  mouth 
of  siren  #2  for  all  conditions. 


80 


Distance  Distance  Distance 


Fig.  A-l.  Sound  Pressure-Time  Function  at  1  kHz  at  Wright- Patterson  Air  Force 
Base  Sound  Facility  in  Anechoic  Mode. 


APPENDIX  B 


HIGH  INTENSITY  TEST  FACILITY  AT  RIVERBANK  ACOUSTICAL  LABORATORIES 

Of  the  three  high  intensity  acoustic  tubes ^^0  available 
at  the  Riverbank  Acoustical  Laboratories,  only  the  low  frequency 
tube  had  adequate  length  to  show  appreciable  nonlinear  propaga¬ 
tion  effects.  Though  originally  designed  for  the  20  to  300  Hz 
region,  operation  to  2000  Hz  was  quite  successful. 

A  drawing  of  the  low  frequency  impedance  tube  is  shown 
in  Fig.  B- 1 .  This  tube  was  mounted  In  the  center  of  a  plywood 
box  filled  with  sand  approximately  8.5  by  12  in.  in  cross  sec¬ 
tion  in  order  to  damp  resonant  vibrations  of  the  tube  walls. 

The  horn  sections  coupling  the  low  frequency  tube  to  the  drivers 
was  also  damped  by  sand.  The  minimum  thickness  of  sand  at  all 
surfaces  of  the  tubes  and  horn  was  4  in.  and  over  6000  lbs.  of 
sand  were  used. 

The  length  of  the  low  frequency  impedance  tube  was  set 
by  the  space  available  for  the  tube  and  its  associated  horn. 

The  tube  was  20  ft.  long,  the  horn  with  its  coupling  section  was 
16.5  ft.  long,  and  the  sample  holder  was  6  ft.  4  in.  long,  giv¬ 
ing  an  overall  length  o£  about  43  ft.  In  the  20  ft.  tube  length 
two  minima  in  the  standing  wave  pattern  could  be  measured  for 
frequencies  above  about  57  cps  for  any  sample  impedance.  Below 
this  frequency,  measurements  were  in  some  cases  restricted  to 
one  minimum  point  and  one  maximum  with  some  loss  in  accuracy  in¬ 
volved  in  calculation  of  the  wave  length  and  estimating  the 
attenuation  along  the  tube.  The  latter  source  of  error  was  not 
very  significant  since  the  attenuation  along  the  tube  was  small 
for  the  samples  of  high  absorption  which  were  measured. 

The  sound  field  inside  the  impedance  tube  was  explored 
by  a  microphone  mounted  inside  the  impedance  tube  at  the  end  of 
a  probe  tube.  A  tube  was  used  instead  of  a  rod  in  order  to  provide 


82 


Impedance  tube  Probe  tub# 


83 


Fig,  B-l.  Low  Frequency  Impedance  Tube. 


a  shielded  conduit  for  the  microphone  cable.  The  probe  tube 
passed  through  a  felt  lined  bushing  at  the  sound- source  end  of 
the  impedance  tube.  The  microphone  end  of  the  probe  tube  was 
supported  by  a  three  legged  carriage  with  felt  at  the  contacts 
between  the  carriage  and  the  inside  surface  of  the  impedance 
tube.  The  felt  at  the  microphone  carriage  supports  and  at  the 
bushing  provided  vibration  isolation  for  the  microphone  which 
is  essential  for  accurate  measurements  of  sound  pressure.  The 
microphones  were  supported  by  their  electrical  connections  (thin 
wires)  to  the  calibrate  resistances  and  to  the  signal  and  cali¬ 
bration  cables  which  were  inside  the  probe  tubes.  The  other  end 
of  the  probe  tube  was  clamped  to  a  movable  carriage  which  sup¬ 
ported  the  microphone  preamplifier. 

The  probe  tube  was  0.75  in.  in  diameter  and  about  23  ft. 
long.  The  preamplifier  carriage  was  moved  along  a  2  by  5  in. 
aluminum  channel  by  means  of  an  endless  chain  driven  by  a  sprocket. 
A  steel  tape  loop  attached  to  the  carriage  was  used  to  determine 
the  microphone  positions. 

The  drivers  for  the  low  frequency  tube  consisted  of  49 
8  in.  cone  loudspeakers,  Jensen  P8-PC.  The  drivers  were  connected 
in  a  seven  by  seven  array  and  each  line  of  seven  drivers  in 
series  was  fused.  It  was  "found  that  the  resonant  frequency  of 
all  the  speakers  was  not  the  same  and,  for  frequencies  near  50  cps, 
high  intensities  could  not  be  attained  without  damage  to  the  cone 
supports.  Under  this  condition,  a  few  speakers  which  were  driven 
near  resonance  had  excessive  amplitudes  and  the  others,  off  reson¬ 
ance,  had  much  smaller  amplitudes.  A  level  of  140  db,  however, 
could  be  obtained  even  in  this  resonance  range. 

The  speakers  were  rated  at  15  watts  for  broadcast  pro¬ 
gram  material.  Using  a  "rule  of  thumb"  conversion  factor  of 
one-half  for  continuous  sinusoidal  signal  would  give  a  rating 
of  7.5  watts  per  speaker.  For  the  nominal  impedance  of  8  ohms, 
the  maximum  current  rating  is  0.96  amp.  For  the  7  by  7  array 
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the  maximum  input  current  is  then  about  7  amps  and  the  maximum 
input  power  is  370  watts.  Operation  at  6  amps  corresponding  to 
290  watts  resulted  in  a  voice  coil  temperature  which  caused  dis¬ 
coloration  of  the  voice  coil  forms,  but  few  failures  occurred 
except  as  noted  for  operation  near  resonance. 

A  pressure  level  of  160  db  corresponds  to  an  intensity 

of  about  one  watt /cm  assuming  that  the  impedance  tube  is  term- 

2 

inated  by  pc.  The  3  in.  diameter  tube  had  an  area  of  45.6  cm 
which  gives  an  acoustic  power  requirement  of  about  45  watts. 

The  efficiency  of  the  driver  and  horn  system  (electrical  input 
power  divided  by  acoustic  power  in  the  impedance  tube)  was  thus 
about  45/290  or  16  percent  over  most  of  the  frequency  range. 

The  speakers  were  mounted  in  a  71  by  69  in.  plywood 
panel  one  inch  thick  and  were  spaced  8.5  in.  on  center.  The 
panel  was  attached  to  the  large  end  of  the  exponential  horn 
which  was  coupled  to  the  low  frequency  tube.  The  horn  construc¬ 
tion  is  described  below.  The  space  behind  the  speakers  was  a 
cavity  of  rectangular  cross  section,  68  by  66  by  8  in.  This 
space  was  ventilated  by  means  of  a  small  blower.  Lined  ducts 
at  the  inlet  and  outlet  openings  for  the  cooling  air  were  used 
to  reduce  the  sound  from  the  back  of  the  speakers. 

The  cone  loudspeakers  were  coupled  to  the  low  frequency 
impedance  tube  by  an  inverted  exponential  horn  designed  to  pro¬ 
vide  an  approximation  to  pc  loading  for  the  loudspeakers.  For 
design  purposes  the  layer  system  was  assumed  to  have  high  absorp¬ 
tion  and  thus  the  impedance  tube  was  terminated  by  a  specific 
impedance  of  about  pc  and  the  specific  impedance  looking  into 
the  source  end  was  also  approximately  the  same.  The  horn  is 
considered  as  a  transformer  which  matches  the  acoustical  imped¬ 
ance  of  the  small  area  at  the  impedance  tube  to  the  large  area 
at  the  loudspeakers.  The  acoustical  impedance  is  defined  as  the 
specific  impedance  divided  by  the  area.  Thus,  a  well  designed 
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horn  will  transform  the  acoustical  impedance  at  the  throat 
pc/S^  to  pc/S2  at  the  large  end  and  give  a  desirable  loading 
for  sound  sources  at  this  end. 

The  horn  was  designed  on  the  basis  of  the  following 
specifications,  a  3  in.  diameter  at  the  throat,  a  59  by  59  in. 
area  at  the  large  end  for  mounting  49  8  in.  loudspeakers,  a 
cut-off  frequency  below  50  cps ,  and  available  space  and  cost 
requirements.  The  equation  for  an  exponential  horn  of  axial 
length  d  is 

exp  (qd)  =  52/8^ 

The  flare  constant  q  determines  the  cut-off  frequency  f^  which 
is  defined  as 

f^  *  qc/47T 

For  d  equal  to  13.6  ft.  and  S2/S^  equal  to  (59/3)^  or  388,  q  is 
0.44  ft" ^  giving  a  cut-off  frequency  of  40  cps. 

In  building  the  horn  it  was  decided,  because  of  cost,  to 
use  flat  plywood  one  inch  thick  for  the  greater  part  of  the  sides 
with  some  deviation  from  an  exponential  expansion  in  area.  The 
length  was  divided  into  seven  segments.  The  first  was  19  in. 
long  and  was  made  of  1/4  in.  sheet  steel  which  was  formed  to  have 
a  circular  cross  section  3  in.  in  diameter  at  the  impedance  tube 
end  and  a  square  cross  section  at  the  other  end.  The  following 
segments  were  made  of  plywood  and  were  each  two  feet  long.  They 
were  hollow  pyramids  with  a  square  cross  section.  At  each  junc¬ 
tion  between  segments,  the  width  of  the  segment  was  fitted  to 
the  exponential  curve  specified  in  the  preceding  paragraph.  The 
large  end  of  the  horn  was  closed  by  a  59  by  59  in.  flat  plywood 
panel  on  which  the  49  loudspeakers  were  mounted  in  a  seven  by 
seven  array  approximately  8.5  in.  on  center . 

The  axis  of  the  horn  was  inclined  from  the  impedance  tube 
axis  as  shown  in  Fig.  B-l.  The  angle  between  the  axes  was  about 
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14°.  A  flanged  steel  tube,  3  in.  I.D.  and  about  3  ft.  long  was 
bent  to  connect  the  impedance  tube  and  horn  so  that  there  was 
no  interference  between  the  microphone  probe  tube  and  its  pre¬ 
amplifier  carriage  and  the  horn. 

Although  the  stepped  horn  differend  from  a  theoretical 
exponential  horn  in  several  respects  and  there  was  an  abrupt 
change  from  the  59  by  59  in.  area  of  the  horn  to  the  smaller 
area  of  the  49  loudspeaker  cones ,  the  horn  did  not  introduce 
pronounced  peaks  and  dips  in  the  frequency  response  of  the  sys¬ 
tem.  The  loading  of  the  loudspeaker  drivers  was  adequate  to 
give  a  fairly  high  efficiency  for  the  conversion  of  electrical 
power  to  acoustic  power  in  the  impedance  tube. 

All  joints  between  the  plywood  sides  of  the  horn  were 
closely  fitted  and  were  glued  and  screwed  together.  The  various 
segments  were  mounted  in  an  enclosure  filled  with  sand  as  des¬ 
cribed  previously. 

The  electronic  equipment  for  the  driving  system  is  listed 
in  Table  B-l,  analyzing  equipment  is  given  in  Table  B-2. 
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TABLE  B-l 


LIST  OF  DRIVING  EQUIPMENT 

Noise  Generator 

Grason-Stadler  Co.  Model  455- B 

Attenuators 

Daven  Attenuators 

0  -  l  db  in  0.1  db  steps,  Type  2506 
0  -  10  db  in  1.0  db  steps.  Type  2511 
0  -  100  db  in  10.0  db  steps,  Type  2513 

Power  Amplifier 

Savage  1  kilowatt  audio  amplifier 
Filter 

General  Radio-- Sound  and  Vibration  Type  1554- E 
Analyzer 
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TABLE  B-2 


LIST  OF  MEASURING  EQUIPMENT 

Low  Frequency  Impedance  Tube  Microphone  Systems 

Massa  M-10I  microphone 

Massa  M- 103  preamplifier 

Massa  M- 104  amplifier  and  power  supply 

Medium  Frequency  Impedance  Tube  Microphone  Systems 

Massa  M-213  microphone 

Massa  M-114  preamplifier 

Massa  M- 185  amplifier  and  power  supply 

Western  Electro-Acoustic  Lab.  Type  100E  amplifier 

Attenuators 

Daven  Attenuators 

0  -  1  db  in  0.1  db  steps,  Type  2506 

0  -  10  db  in  1.0  db  steps,  Type  2511 
0  -  100  db  in  10.0  db  steps,  Type  2513 

Filters 

Allison  variable  low  pass  and  high  pass  filters 
Analyzer 

Quan-Tech  Model  303 
Voltmeter 

Ballantine  electronic  voltmeter,  Model  300 


89 


Unclassified _ 

Security  Classification 


DOCUMENT  CONTROL  DATA  •  R&D 

(Security  claaaittcation  ot  title,  body  ot  abatratt  and  mdantng  animation  mm*  bp  entered  *h»n  the  overall  import  ta  claaeit(ed) 


1  ORIGINATING  ACTf^l^Y  (Corporate  author) 

Za  REPORT  IECUWITY  CLASSIFICATION 

IIT  Research  Institute 

Unclassified 

10  West  35  Street 

Zb  GROUP 

Chicago,  Illinois  606l6 

J  REPORT  TITLE 


Propagation  of  High  Intensity  Sound  Waves 


4  DESCRIPTIVE  NOTES  (Type  o(  report  and  fncJuaft'e  dates) 

Final  January  1966  to  March  1967 

T” AUTHORS  (Lamt  name,  tint  name.  Initial) 

Karplus,  Henry  B. ,  Raelson,  Vemer  J.,  Hruska,  Gale  R. 


G.  REPORT  DATE 

September  19 67 

?  a  TOTAL  NO.  OF  PAGE* 

89 

7b.  NO-  OF  R EFI 

24 

G*  CONTRACT  OR  GRANT  NO. 

9  a  ORIGIN  A  TOR* 1*  REPORT  NUMBERS.) 

AF  33(6l5)-3436 

affdl-tr-67-37 

t>,  PROJECT  NO 

4437 

C- 

Task  443701 


0  b.  other  REPORT  NCtfSj  (Any  other  number*  (hat  may  ba  aaeigned 
thte  report} 


I  a-  AVAIL  A0ILITY/LIHIT  A  TIQN  NOTICES™*  Hnf'nTnpnt 

Distribution  limited  to  U.  S.  GovemmeSf  ™  , 

agencies  only;  test  and  evaluation; 

^statement  applied  r** /+  ?/ 

1  Other  requests  for  this  do^11  snt  must  bo 
referred  to  AF  Flight  D^n^i***-*  '.abora.ory, 

Wright-Patferson  AF6*  Ohio  45433 


ia,  sH*E>n  soring  military  activity 

Air  Force  Flight  dynamics  Laboratory 
Wright-Patterson  AFB,  Ohio  45433 


13  ABSTRACT 

Sound  field  characteristics  were  studied  as  a  function  of  increasing 
intensity.  It  was  determined  that  one  effect  out  weights  others.  This  is  the 
distortion  generated  within  the  siren  which  makes  it  impossible  to  generate 
pure  tones  above  a  certain  intensity.  At  high  intensity  the  waveform  of  the 
sound  assumes  a  sawtooth  shape  producing  periodic  shock  waves.  Conversely ,  it 
is  concluded  that  if  the  distortion  is  small  at  the  mouth  of  the  siren,  additional 
distortion  in  the  subsequent  propagation  will  not  be  significant,  Waveforems 
observed  in  the  AFFDL  facility  are  included  in  an  appendix. 


Interaction  of  sound  waves  from  different  sirens  will,  in  addition,  give 
rise  to  intermodulation  effects  producing  sum  and  difference  frequencies.  In 
general,  the  intensity  in  the  region  in  which  the  radiation  from  different  sirens 
can  interact  will  be  so  much  lower  than  the  intensity  in  the  throat  of  the  siren 
horn  that  these  effects  will  also  be  small  compared  with  harmonic  distortion 
generated  in  the  sirens. 


DISTRIBUTION  OF  THIS  ABSTRACT  IS  UNLIMITED 


DD  1473  _ Unclassifi ed _ 

Security  Classification 


Unclassified 

Security  Classification 


LINK  a 


KEY  WORDS 


Acoustic  Environment 
Sound  Wave  Propagation 


INSTRUCTIONS 


1*  ORIGINATING  ACTIVITY:  Enter  the  name  an d  address 
of  the  contractor,  subcontractor,  grantee,  Department  of  De¬ 
fense  activity  or  other  organization  ( corporate  author)  issuing 
the  report. 

2a.  REPORT  SECURITY  CLASSIFIC ATlON:  Enter  the  over¬ 
all  security  classification  of  the  report.  Indicate  whether 
“Restricted  Data*1  is  included.  Marking  is  to  be  in  accord¬ 
ance  with  appropriate  security  regulations, 

2 h.  GROUP;  Automatic  downgrading  is  specified  in  DoD  Di- 
rective  5200.  10  and  Armed  Forces  Industrial  Manual.  Enter 
the  group  number.  Also,  when  applicable,  show  that  optional 
markings  have  been  used  for  Group  3  and  Group  4  as  author¬ 
ized. 

3*  REPORT  TITLE:  Enter  the  complete  report  title  in  all 
capital  letters.  Titles  in  all  cases  should  he  unclassified. 

If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  show  title  classification  in  all  capitals  in  parenthesis 
immediately  following  the  title, 

4.  DESCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e . g. ,  interim,  progress,  summary,  annual,  or  final. 

Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHOR(S):  Enter  the  natne(s)  of  authors)  as  shown  on 
or  in  the  report.  Enter  last  name,  first  name,  middle  initial. 

If  military,  show  rank  branch  of  service.  The  name  of 
the  principal  »:  t hor  is  an  absolute  minimum  requirement. 

6.  REPORT  DATL:  Enter  the  date  of  the  report  as  day, 
month,  year,  or  month,  year.  If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7a.  TOTAL  NUMBER  OF  PAGES:  The  total  page  count 
should  follow  normal  pagination  procedures,  L e. ,  enter  the 
number  of  pages  containing  information. 

7ft.  NUMBER  OF  REFERENCES:  Enter  the  total  number  of 
references  cited  in  the  report, 

Sa.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written, 

Sb,  8c,  At  fid.  PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification*  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9a.  ORIGINATOR'S  REPORT  NUMBER(S):  Enter  the  offi¬ 
cial  report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.  This  number  must 
be  unique  to  this  report, 

9b.  OTHER  REPORT  NUMBER(S):  If  the  report  has  been 
assigned  any  other  report  numbers  (either  by  the  originator 
or  by  Che  sponsor),  also  enter  this  number(s). 

10.  .'VAIL ABILITY/LIMITaTION  NOTICES:  Enter  any  lii 
itations  on  further  dissemination  of  the  report,  other  than  thf-se 


imposed  by  security  classification,  using  standard  statements 
such  as: 

(!)  “Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC," 

(2)  “ Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized," 

(3)  “U.  S,  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC,  Other  qualified  DDC 
users  shall  request  through 


(4)  “U.  S.  military  agencies  may  obtain  copies  of  this 

report  directly  from  DDC.  Other  qualified  users 
shall  request  through 

it 


(5)  "All  distribution  of  this  report  is  controlled.  Qua! 
ified  DDC  users  shall  request  through 

it 


If  the  report  has  been  furnished  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  fact  and  enter  the  price,  if  known, 

H.  SUPPLEMENTARY  NOTES:  Use  for  addition^  explana¬ 
tory  notes. 

12.  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  (pay¬ 
ing  (or)  the  research  and  development.  Include  address. 

13.  ABSTRACT:  Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  the  body  of  the  technical  re¬ 
port.  If  additional  space  is  required,  a  continuation  sheet  shall 
be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  reports 
be  unclassified.  Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  military  security  classification  of  the  in¬ 
formation  in  the  paragraph,  represented  as  (TS)r  ( S),  fC>.  or  (U j 

There  is  no  limitation  on  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  150  to  225  words 

14.  KEY  WORDS:  Key  words  are  technically  meaningful  terms 
or  short  phrases  I  hat  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.  Key  words  must  he 
selected  so  that  no  security  classification  is  required.  Identi¬ 
fiers,  such  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  he  followed  by  an  indication  of  technical  con¬ 
text.  The  assignment  of  links,  rules,  and  weights  is  optional. 


Unclassified 


Security  Classification 


